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ABSTRACT 


Several  milestones  ha\-e  been  achieved  tou-ards  the  realization  of  a  practical  space  plasma 
thruster  from  the  tandem  mirror  rocket  experiment,  i.e.  the  specifi  impulse,  thrust,  energj' 
conversion  efficiency,  and  mass  flow- rate  have  been  determined.  The  experiment  operates  at 
9.4  k\V  of  input  power  at  an  rf-to- plasma  efficiency  of  68%;  The  ion  temperature  is  172  e\' 
(  2,000,000  K  )  which  gives  an  I^p  of  12,852  s.  The  thrust  is  76  mN  (milli  New-ton)  which  is 
a  high  \’alue  for  a  low  input  power  (9.4  kW)  and  very  high  (  12,852  s  ).  The  radiation 
loss  was  found  to  be  very  low.  Most  important,  these  results  fall  within  our  prediction. 

A  key  advantage  of  this  propulsion  system  is  its  high  energy  density  that  can  be  contained 
in  the  plasma.  A  flight  system  design  and  analysis  has  shown  that  the  present  experimental 
device  can  handle  a  power  level  in  excess  of  10  MW,  resulting  in  specific  weights  as  low  as  0.04 
kg/kW.  The  overall  specific  weight  of  a  space  vehicle  for  a  piloted  Mars  mission,  including 
nuclear  power  and  power  conditioning  equipment,  is  only  8  kg/kW.  This  is  a  conser\’ative 
value  which  can  be  further  reduced  to  6  kg/kW. 

A  1 50  day  trip  to  .Mars  (  180  day  round  trip  including  90  days  stay  over  )  ran  be  achieved 
with  <.nly  10  MW  of  power.  The  trip  time  has  been  optimized  by  varying  the  and  thrust 
to  meet  a  required  returned  mass  of  .57.000  kg.  In  that  study  an  efficiency  of  only  50%  was 
assumed.  Except  for  nuclear  power,  the  technologj-  involved  is  state-of-the-art  with  very 
minimal  engineering  scaling.  In  addition,  substantial  ad\'ancements  in  high  temperature 
superconductors,  solid  state  transmitters  and  pow’er  conditioning  equipment  are  expected. 
Such  development  makes  this  device  a  most  competitive  electric  propulsion  approach  when 
it  is  fully  developed.  The  concept  can  be  demonstrated  near-term  as  gauged  on  our  rapid, 
limited  resource  progress. 

Our  future  goal  is  to  reach  a  plasma  coupling  efficiency  as  high  as  90%.  This  expectation 
is  reasonable  because  there  is  sufficient  room  for  improvement  as  shown  by  our  theoretical 
emalyses.  These  studies  predict  that  higher  coupling  efficiency  can  be  obteiined  if  the  mag¬ 
netic  field  is  made  mure  unifurin  and  the  plasma  radius  larger.  At  present,  the  ionization 
efficiency  is  very  low  because  the  microwave  frequency  is  not  in  resonance  with  electron  cy¬ 
clotron  frequency.  The  efficiency  is  expected  to  improve  when  the  newly  installed  14.5  GHz 
microwave  transmitter  becomes  operational.  This  frequency  is  closer  to  resonance  at  the 
operating  fields  of  the  experiment  (  2.0  to  10.0  kG  ).  These  features  will  be  implemented  in 
our  future  experiment. 

A  generaJ  review  of  this  research  is  given  in  the  introduction  section.  Section  2  summa¬ 
rizes  the  progress  made  in  the  past  year.  The  detailed  discussions  of  the  theoreticail  work, 
thrust  balance  development,  laser  fluorescence  system  development,  flight  system  design,  ex¬ 
haust  modeling  and  experimental  results  are  given  in  sections  3,  4,  5,  6  and  7.  Among  these 
topics,  the  theoretical  analysis  of  wave  propagation,  some  important  experimental  results 
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and  flight  sv>item  dp'^ign  and  exhatist  modeling  have  been  presented  in  the  previoiis  reports, 
riiev  are  in*lu<!t  d  here  to  give  the  complete  picture  and  to  show  the  depth  and  I  '• 
of  this  research  project. 


1  Introduction 

We  have  often  pointed  out  that  the  tandem  mirror  plasma  rocket  has  the  potential  of  being  a 
highly  efficient  propulsion  device  with  vTiriable  and  thrust  from  low  to  very  high  values. 
Such  a  propulsion  system  hinges  on  the  reliable  production  and  controlled  exhaust  of  plasma 
from  an  rf  heated  tandem  mirror  device.  To  this  general  end,  an  experimental  research 
program  was  initiated  at  MIT  in  1988  [1],[2],  [3], (4], '5], [6]. 

It  is  the  point  of  the  time  and  for  future  development  that  we  reiterate  the  advantages  of 
this  rocket  engine,  as  weel  as  the  reasons,  purposes  and  scenario  of  this  research.  It  has  been 
shown  in  fusion  energy  research  on  a  Tandem  Mirror  device  that  the  plasma  can  be  heated  to 
as  high  3,000.000  K,  but  at  very  low  density  ( typically  in  the  neighborhood  of  5  x  10^*m“^). 
This  high  plasma  ion  temperature  will  give  a  very  high  specific  impulse  not  attainable  by 
other  electrical  propulsion  methods  with  the  exception  of  ion  propulsion.  The  thrust  of  an 
ion  engine  is  very  low  while  the  thrust  of  a  tandem  mirror  engine  can  be  potentially  very 
high  by  operating  the  system  at  low  temperature  and  high  density.  A  very  important  feature 
unique  to  this  rocket  engine  is  that  both  the  thrust  and  specific  impulse  can  be  v-aried  by 
controlling  the  temperature  and  density.  Also  the  relationship  of  the  variation  of  these  two 
parameters  in  this  engine  is  the  most  desirable,  i.e..  this  engine  can  be  operated  at  high 
thrust  and  low  specific  impulse  and  vice  versa  at  constant  power  level,  whereas  in  the  other 
types  of  electric  systems  high  thrust  always  associates  with  high  specific  impulse  and  high 
power.  Therefore  high  thrust  ran  be  delivered  to  the  vehicle  when  it  is  accelerated  from  low 
speed  and  the  thrust  can  be  tapered  off  when  the  vehicle  is  approaching  a  cruising  speed.  The 
spacecraft  can  then  be  flown  at  constant  power.  This  is  an  economical  way  of  utilizing  the 
precious  space  power  auid  an  efficient  way  of  flying.  Another  advantage  is  that  this  system 
can  be  operated  at  either  steady  state  or  pulse  mode.  The  energy  conversion  efficiency  is 
very  high  with  ICRF  (ion  cyclotron  resonance  frequency  )  heating  method. 

However,  the  tandem  mirror  device  is  by  no  means  mature  enough  for  direct  application 
to  space  propulsion.  As  it  was  poi.ited  out  previously,  the  Tandem  Mirror  devices  for  fusion 
energy  research  are  lairge  and  heavy,  where2a  the  propulsion  device  must  be  compact  and 
light  weight.  Hence  we  must  first  demonstrate  that  the  projected  plasma  parameters  are 
achievable  in  a  compact  device  (  which  was  accomplished  in  the  first  year  of  our  research 
).  For  propulsion  purposes  the  preferred  plasma  temperature  is  much  lower  than  the  tem¬ 
perature  required  for  ignition,  with  much  higher  density.  Low  temperature  and  high  density 
plasma  conditions  have  to  be  established.  The  mode  of  operation  is  also  different.  For  a 
fusion  energy  reactor  the  plasma  has  to  be  well  confined  and  the  end  leakage  should  be  kept 
as  low  as  possible.  On  the  other  hand,  for  propulsion  purposes  the  fuel  gas  has  to  be  fed. 
ionized,  heated  and  exhausted  in  a  controlled  feishion.  Such  a  dynamic  mode  of  streaming 
plasma  operation  is  unknown  and  hzis  to  be  developed. 
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The  devirahle  temperature  range  for  propulsion  application  is  from  5  to  100  e\’  (1.16 
iiiilli'in  (Kgrees  K)  and  density  range  is  fr<im  10'*cm  ^  to  10*“cm“^.  To  obtain  such  a  wide 
range  of  plasma  parameters  extensive  research  is  required.  The  physics  of  wave  propaga¬ 
tion  and  dissipation  in  the  plasma  in  a  mirror  geometry  had  not  yet  been  fully  investigated 
before  the  mirror  fusion  program  was  terminated.  Since  the  ICRF  heating  is  chosen  as  the 
method  for  transfering  electric  energy  to  the  propulsive  energy,  a  full  understanding  of  the 
wave-plasma  physics  has  to  be  realized.  The  management  of  the  advanced  propulsion  group 
at  JPL  appreciated  this  need  and  wisely  advised  us  to  concentrate  on  the  understanding  of 
the  physics  involving  the  ICRF  heating  at  the  time  when  this  project  was  funded.  Therefore 
analytic,  numerical  and  experimental  studies  on  the  physics  of  wave  propagation  and  ab¬ 
sorption  in  the  plasma  of  the  mirror  device  have  been  undertaken  as  our  second  major  task. 
This  effort  has  been  paid  off  and  a  sufficient  understanding  of  the  physics  of  ICRH  has  been 
obtained  and  is  the  key  to  the  success  of  our  experiment.  This  leads  to  the  achievement  of 
stable  plasma  condition  and  heating  of  ion  temperature  to  170  keV'. 

To  make  this  system  applicable  to  space  propulsion  the  plasma  has  to  be  detached  from 
the  field  lines  after  being  exhausted  through  one  end  while  the  other  end  is  closed  or  used 
for  fueling.  The  dynamic  property  of  the  plasma  for  the  open  ended  operation  with  fueling 
and  exhaust  are  unknown.  Two  plasma  detachment  schemes  have  been  proposed;  One  is  the 
disconnection  of  magnetic  field  line  ,  and  the  other  is  the  injection  of  annular  cold  geis  (  a 
hot  plsisma  exhaust  surrounded  by  a  layer  of  cold  gas  named  as  hybrid  plume  ).  Theoretical 
in\estigation  on  the  understanding  of  the  dynamic  property  of  plasma  exhaust  has  been  an 
<»n  going  effort  since  conceiving  this  concept  in  early  1980.  The  detachment  of  the  plasnja 
with  the  annual  injection  of  gas  has  been  found  feasible.  .At  present  the  experiment  can  not 
been  <lone  because  the  fueling  and  exhaust  systems  ha\'e  not  been  implemented  due  to  lack 
of  funding.  However,  the  thrust  and  mciss  flow-rate  have  been  measured  at  the  exit  of  the 
exhaust  chamber  and  found  to  be  in  agreement  with  prediction.  Therefore  the  next  major 
task  will  the  experimental  demonstration  of  the  plasma  detachment. 

Although  the  spacecraft  system  design  and  the  analysis  of  engine  perfomance  in  a  model 
mission  have  never  been  part  of  our  research  program,  we  are  fully  aware  of  the  importance 
of  the  system  implementation.  Therefore  a  scoping  study  of  the  system  design  and  mission 
analysis  based  a  assumed  trip  to  Mars  has  been  carried  out.  From  this  exercise  we  learned 
that  this  system  is  very  competitive  in  terms  of  specific  weight  and  trip  time. 

This  device  is  very  compact.  It  is  only  three  meters  in  length  and  less  than  one  meter  in 
diameter.  It  is  only  one  tenth  of  the  size  of  the  smallest  fusion  device  [7]  but  at  four  time 
the  field  intensity.  The  machine  parameters  are  shown  in  Table  1  and  the  system  layout  is 
shown  in  Figure  1.  Although  the  device  is  much  smaller  than  the  fusion  device,  the  needed 
experimental  facility  is  as  extensive  :ts  that  for  fusion  research.  The  present  size  w-as  chosen 
based  on  the  minimum  spaces  needed  for  diagnostic  purpose.  In  reality,  the  actual  size  can 
be  much  smaller.  The  system  scoping  study  indicates  that  this  device  can  deliver  10  MW  of 
power  which  is  large  enough  for  a  Mars  mission. 
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Tabic  1:  Device  Paraniclcrs 


Overall  length  |  102  m 

Central  cell  length  :  1.15m 
Central  cell  radius  i  0.36  rn 
Central  cell  field  (Be)  j  2.4  kG 
Maximum  field  (Bmax)  |  20.0  kG 
Plasma  heating  power:  j 

•  ECRF  j  1.0  kW 

•  ICRF _ I  1.0  MW 


The  experiments  conducted  on  this  device  have  shed  new  light  on  the  dynamics  of  these 
high-power  devices  and  their  application  to  multimegawatt  space  propulsion.  In  the  past 
years  the  successful  ICRF  heating  of  a  mirror  plasma  has  been  demonstrated  both  theoreti¬ 
cally  and  experimentally  ;8],  91.  Two  types,  I  and  II,  of  plasma  discharges  can  be  produced, 
with  greater  stability  being  obtained  in  type  II  at  high  magnetic  fields.  Plaisma  parzuneters 
nearly  equiA-alent  to  those  of  large  fusion  devices  have  been  obtained. 

Further  milestone  accomplishments  are  the  measurement  of  ion  temperature,  the  mea- 
.surement  of  thrust  and  mass  flow-rate.  The  ion  temperature  is  an  important  pzirameter 
because  it  determines  b<ith  and  thrust.  The  ion  temperature  has  not  been  measured  in 
any  electric  propulsion  system.  The  radiation  loss  was  found  to  be  very  low.  The  extracted 
efficiency  was  as  high  as  68*^  even  though  the  rf  transmitter  was  running  at  a  deteriorated 
condition.  At  an  /,p  of  12,852  s  and  p<»wer  level  of  t>nly  9.4  k\V  the  thrust  of  76  mN  is  a 
respectable  value.  Most  important,  these  results  fall  within  our  prediction.  These  achieve¬ 
ments  are  way  ahead  of  the  goals  which  we  set  out  to  accomplish. 

The  10  MW  class  vehicle  and  MARS  mission  were  chosen  as  the  models  for  our  flight 
system  and  mission  analysis  so  that  a  comparison  can  be  made  with  an  established  anzdysis 
for  an  ion  engine.  The  design  study  of  a  20  kW  class  Tandem  Mirror  Rocket  based  on  the 
fuel  cell  or  solar  panel  technology  is  undervray.  A  direct  comparison  with  MPD  and  Arcjet 
can  then  be  made  when  such  a  system  design  is  ready.  The  10  MW  class  system  study 
is  comparable  to  the  ion  engine  and  is  expected  to  show  better  performance  w’hen  higher 
efficiency  and  lower  weight  are  considered.  In  addition  to  high  /,p  and  efficiency,  this  rocket 
has  the  flexibility  of  being  operated  at  either  variable  or  constant  and  thrust.  There  is 
no  electrode  damage  problem  and  the  rocket  can  use  any  type  of  propellant. 

Our  future  goal  is  to  reach  an  efficiency  as  high  as  90%.  This  expectation  is  reasonable 
because  there  is  sufficient  room  for  improvement  as  shown  by  our  theoretical  analysis.  The 
newly  accomplished  theoretical  analysis  predicts  that  higher  coupling  efficiency  can  be  ob¬ 
tained  if  the  magnetic  field  is  made  more  uniform  and  the  plasma  radius  larger.  At  present, 
the  ionization  efficiency  is  very  low  because  the  microwave  frequency  is  not  in  resonance  wuth 
electron  cyclotron  frequency.  The  efficiency  is  expected  to  improve  when  the  newly  installed 
14.5  GHz  microwave  transmitter  becomes  operational.  These  features  will  be  implemented 
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Figure  1.  Experimental  set 


in  our  ftiture  experiment  and  will  undoubtedly  raise  the  efficiency.  The  weight  for  the  power 
Conditioning  ecjuiptnent  was  scaled  linearly  fr*(jn  a  1  MW  system.  Actually  the  specific 
weight  will  level  off  when  the  power  increases.  The  frequency  of  the  transmitter  system  used 
as  basis  for  the  system  study  is  100  MHz  which  is  much  higher  than  the  500  kHz  to  10  MHz 
frequency  needed  for  propulsion  application.  The  weight  of  the  power  conditioning  system 
will  be  much  lighter  at  these  frequencies.  The  projected  specific  weight  can  be  reduced  t<f 
as  low  as  6  kg/kW.  These  values  are  significantly  lower  than  any  known  electric  propulsi<^n 
system. 

To  bring  our  system  into  the  mainstream  of  electric  propulsion  programs,  it  is  our  wish  to 
to  move  forward  to  carry  out  the  experiment  that  will  demonstrate  the  plasma  detachment. 
However,  such  an  experiment  requires  a  major  upgrade  of  our  facility  which  is  not  realistic 
at  the  current  funding  level.  We  propose,  instead,  to  carry  out  the  experiments  to  improve 
the  energy  conversion  efficiency  and  obtain  the  conditions  for  high  thrust  operation.  The 
development  of  reliable  plasma  condition  control  is  as  important  as  the  detachment.  The 
detachment  may  become  obvious  at  higher  thrust  (  high  density  and  low  temperature  ). 

The  entire  facility  hcis  been  built  and  operated  by  graduate  students,  undergraduate 
students,  part  time  staffs,  and  visiting  scientists.  Of  the  the  entire  facility  only  the  mirror 
device  was  fabricated  with  the  instrumentation  fund  from  the  Air  Force  and  constitutes  a  very 
small  fraction.  The  major  components  of  the  facility  are  surplus  equipment.  This  project  has 
made  a  very  significant  contribution  to  education.  There  are  two  graduate  students  earned 
their  doctoral  degrees  based  on  their  work  on  this  project.  In  addition  there  have  been 
four  graduate  students,  five  undergraduate  students  and  even  summer  high  school  students 
working  on  this  project.  This  project  has  been  greatly  benefited  from  the  visiting  scientists, 
two  of  them  from  the  Institute  of  Physics  of  the  Chinese  Academy  of  Science  in  Beiging. 
one  from  the  Southwestern  Institute  of  Physics  in  Sichuan,  China.  They  have  helped  in  the 
constriiction  of  the  facility,  developed  the  laser  fluorescence  system,  many  probes,  and  the 
development  of  space  bound  solid  state  RF  transmitter  systems. 
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2  Summary  of  Progress 

The  progress  made  in  the  past  year  is  summarized  in  the  following  subsections.  To  keep  the 
discussion  brief,  only  the  key  points  will  be  given  and  details  can  be  found  in  later  sections. 

2.1  Experimental  Study 

The  experimental  setup  is  shown  in  Figure  2.  A  wave  is  launched  from  a  double-half- loop 
antenna  located  at  the  high  field  region  of  the  central  cell.  The  wave  propagates  radially 
inwards  and  toward  the  resonance  midplane  where  the  diagnostics  are  located.  A  172  e\' 
ion  temperature,  a  21.6  eV  electron  temperature  and  a  1.98  x  10^^m“^  plasma  density  are 
obtained  from  probe  measurements.  In  addition,  the  absence  of  pronounced  Ha  lines  in  the 
Doppler  broadening  spectrometer  confirms  the  expectation  of  a  high  degree  of  ionization  and 
hence  a  relatively  low  power  loss  due  to  atomic  excitation.  Other  important  parameters  can 
be  extracted  from  these  data  and  are  included  in  Table  2.  This  operating  point  is  shown 
in  the  performance  graph  of  Figure  3  at  the  power  level  of  9.4  kW.  From  these  data,  the 
rf-to-plasma  coupling  efficiency  is  determined  to  be  68%. 


Table  2:  Plasma  Parameters 


T, 

Measured  ion  temperature 

172  eV  1 

Te 

Measured  electron  temperature 

21.6  eV  i 

■Measured  plasma  density 

1.98xl0*^m  3 

T 

Extracted  energy  confinement  lime 

50  p  s 

P.n(rf) 

Extracted  input  rf  power 

j  9.4  kW 

/in(microw’ave) 

'  Extracted  input  microwave  power 

i  0.3  kVV  : 

Pab»orb(i*^n) 

Power  absorbed  by  the  ions  Upk  TAp 

5.2  kW 

!  Pa6sorb(electron)  :  Power  absorbed  by  the  electrons  Uek  Telp  j  0.69  kW 

Pexhaust  !  Power  flows  iuto  the  exhaust  chambers  |  0.5  kW 

^  Coupling  efficiency  |  68% 

lap  Extracted  specific  impulse  ,  12,852  s 

m  Mass  flow’- rate  }  0.1  mg  ; 

Thrust  j  Measured  I  76  mN  ' 

Resource  limitations  have  prevented  us  from  operating  in  the  streaming-plasma  mode. 
This  capability  requires  that  a  vacuum  interface  to  the  large  vacuum  tank  in  our  laboratory 
be  manufactured  and  the  design  and  manufacture  of  a  propellant  injector.  Tw’o  additional 
magnet  power  supplies  are  required  to  run  the  system  asymmetrically.  These  system  en¬ 
hancements  are  necessary  to  perform  the  plasma  detachment  experiment. 

Nevertheless,  rapid  progress  has  been  made  in  the  understanding  of  both  the  plasma  and 
neutral  components  in  the  exhaust.  The  tandem  mirror  is  operated  at  a  mirror  ratio  of  10 
and  the  plasma  is  well  confined.  .Any’  exhaust  plasma  is  the  result  of  the  loss-cone  effect 
and  is  only  20%  of  the  central  density.  However,  this  axial  loss  plasma  shows  no  drop  in 


Lnser  Fluorescence 
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Microwave  Interferometer 


1000 

Figure  3:  Overlay  of  the  experimental  data  with  performance  projection 
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temperature  from  that  of  the  central  cell.  Radially,  both  density  and  temperature  profiles 
are  nearly  flat  suggesting  no  loss  of  efficiency  due  to  profile  effects.  A  laser  fliK»rescence 
system  is  now  operating  to  measure  neutral  density.  The  system  utilizes  a  new  two-t>ptical- 
path  fluorescence  signal  reduction  method  developed  by  our  group  10  .  The  radial  profile  <jf 
the  neutral  density  in  the  exhaust  has  been  measured.  The  results  show  a  hollow  profile, 
peaking  near  the  edge  and  dropping  to  negligible  values  on  axis.  This  phenomenon  agrees 
with  our  earlier  prediction  that  the  wall  is  insulated  by  neutrals.  The  density  in  the  exhaust 
will  match  that  of  the  central  cell  when  the  device  is  operated  in  the  streaming-plasma  mode. 

Two  independent  thrust  me«isurement  methods  are  being  developed:  a  grid  energy  an¬ 
alyzer  to  measure  ion  axial  energy  from  which  the  plasma  velocity  and  thrust  distribution 
can  be  determined,  and  a  ballistic  balance  to  measure  total  thrust.  The  former  method 
requires  a  sweeping  power  supply  or  amplifier  which  has  been  acquired  and  will  be  ready  for 
the  future  experiment.  The  ballistic  balance  has  been  successfully  used  to  measure  the  the 
thrust  of  this  device.  The  sensitivity  of  the  balance  is  better  than  0.1  mN.  The  measured 
thrust  at  exhaust  is  given  in  Table  2.  The  description  of  the  ballistic  balance  are  given  in 
section  4. 

The  Doppler  broadening  measuring  instrument  is  a  monochromatic  spectrometer  fitted 
with  a  fast  rotating  mirror  in  the  optical  path.  The  system  was  developed  in  the  past  few 
years  and  has  become  operational  recently.  A  resolution  as  high  as  0.5  .Angstrom  has  been 
achieved.  However,  we  have  difficulty  in  making  measurement  because  the  radiation  was 
found  to  be  extremely  weak.  This  is  a  strong  indication  of  low  radiation  losses  from  the 
plasma. 

2.2  Theoretical  Analysis 

.As  has  been  reported  in  the  past,  the  experimental  effort  was  guided  by  sound  theoretical 
predictions.  The  ezirly  theoretical  work  was  the  development  of  the  numerical  simulation 
BEACH  code  ^8].  In  the  past  year  an  analytical  solution  of  ICRF  wave  propagation  in  an 
axially  inhomogeneous  cylindrical  plasma  has  been  undertaken [9^.  The  main  result  of  this 
work  has  been  the  development  of  the  code  CA'^LWAVE.  This  new  tool  allows  us  to  carry 
out  some  parameter  scans  which  can  not  be  done  with  numerical  simulation.  The  results 
reveal  two  important  points.  First,  the  behavior  of  the  plasma  impedance  as  a  function  of 
the  mirror  ratio  is  in  good  agreement  wdth  the  resonance  volume  theory.  This  theory  sim¬ 
ply  states  that  the  dissipation  is  proportional  to  the  volume  of  plasma  that  is  close  to  the 
resonance.  Second,  optimal  power  coupling  to  the  plasma  is  obtained  with  the  resonance 
located  directly  underneath  the  rf  antenna  or  at  the  machine  midplane.  These  findings  pro¬ 
vide  appropriate  design  guidelines  which  favor  machines  providing  a  large  plasma  volume  in 
a  uniform  magnetic  field  at  the  ion  cyclotron  resonance. 

The  results  obtained  for  the  slab  geometry  show  that  the  efficiency  T  increases  as  reso¬ 
nance  is  approached.  In  other  words,  large  wave  power  is  dissipated  or  the  wave  is  damped. 


These  theoretiral  predictions  agree  with  preliminary  experimental  results  of  the  axial  B, 
profile  as  have  been  reported  previously. 

These  results  will  guide  our  future  experiments  even  though  many  hardware  nK>difications 
will  have  to  be  made  in  order  to  implement  the  new  features. 


2.3  Flight  System  Study 

A  flight  system  design  study  has  been  undertaken  (postulating  a  piloted  Mars  mission  at 
10  -VlWe  and  50%  efficiency )[  11].  The  preliminary  results  are  very  attractive.  The  overall 
specific  weight  of  a  10  MW  flight  system  is  only  8  kg/kW  using  an  existing  space  nuclear 
electric  system  design.  Initial  calculations  performed  recently  for  our  group  by  J.  Gilland 
and  C.  Hack  of  N.ASA/Lewis  |12i  indicate  an  attractive  150  day  transit  to  MARS.  The  total 
trip  time  is  480  days  including  a  90  day  stay  on  MARS.  The  calculations  used  an  input  value 
of  8  kg/kW  and  an  efficiency  of  50%.  Since  then,  experimental  results  and  hardw'are  design 
modifications  have  allowed  us  to  refine  the  subsystem  specific  weight  estimates,  as  well  as  to 
consider  a  more  realistic  efficiency  value  of  70%. 
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3  Theoretical  Analysis 

The  effort  of  the  theorel  ical  work  is  to  obtain  much  more  rigorous  treatment  (»f  the  ICRF  wave 
propagation  than  previously  reported.  The  true  geometry  can  be  considered  in  numerical 
simulations,  but  the  physical  insight  is  difficult  to  extract.  The  functional  dependence  of 
the  wave  dissipation  on  the  plasma  parameters  can  be  formulated  analytically  only  after 
many  simplified  assumptions  are  made.  Notwithstanding  these  simplifications  building  upon 
the  relatively  crude  analytical  solution  of  slab  and  cylindrical  models  which  was  presented 
previously  [6].  The  refined  analysis  of  both  models  was  undertaken  and  the  results  of  the 
slab  model  is  discussed  in  the  following. 


3.1  Analytical  Study 

3.1.1  Slab  Model 


The  analysis  of  wave  propagation  in  the  slab  geometry  is  much  simpler  than  that  in  the 
cylindrical  geometry  and  can  provide  information  on  the  coupling  efficiency  of  ICRF  power 
to  plasma.  The  detailed  analysis  will  be  given  in  reference  [131.  A  figure  of  merit  was  defined 
as  the  dissipated  power  divided  by  the  source  power  as  given  by  the  incident  Poynting  Flux. 
Since  for  the  majority  of  the  cases,  u;  —  fij  3>  the  much  simplified  ^  is  obtained 


^  =  2  cos  dr 


C) 


cos*  9j 

{cosOj  -J-  ^  cos^t)^  +  sin*  0/J 

where  dij  are  the  incident  and  transmitted  angles  of  the  power  flux  and  are  defined  as 


(1) 


9, 

9f 


where 


S^-D^ 

s 

,5*-D* 

We  will  examine  the  following  two  important  limits: 


kla  -  kl 


1.  Resonance  limit  (u;  — >  fij 

In  this  limit,  it  is  seen  that  5,  D  — »  ex;.  This  allows  the  dropping  of  the  last  term  in 
the  denominator.  It  is  also  seen  that  cos  dr  —  1,  and  kpi  ':$>  ko.  This  results  in  the 


7^- 


(2) 


Since  rr  1,  Dfie  r:  0  in  the  limit  of  resonance,  it  can  be  shown  that  the  efficiency 
^  approaches  a  maximum  of  1  when  resonance  is  approached. 
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2.  Far  from  the  Resonance  Limit  (u;  »  Q;,  wi;  <g;  fi,) 


For  frequencies  far  away  from  resonance,  the  value  of  the  efficiency  can  be  given  by  the 
same  equation  as  before  (Equation  2).  However,  the  simplification  of 


f  ^  Sfu 


V  “  ^Re 

can  still  be  made.  This  results  in  an  efficiency  of 


»(2) 


cos  dj 


(3) 


(^) 


Since  for  our  range  of  frequencies,  flj  for  frequencies  off  resonance,  the  efficiency 

approaches  zero. 


The  results  obtained  for  the  slab  geometry  show  that  J-  increases  as  resonance  is  ap¬ 
proached.  In  other  words,  large  w’ave  power  is  dissipated  or  the  wave  is  damped.  These 
theoretical  predictions  agree  with  preliminary  experimental  results  of  the  axial  Bz  profile. 


3.1.2  Cylindrical  Model 

Purpose  of  Study  -  The  purpose  of  this  study  is  to  examine  the  problem  of  ICRF  wave 
propagation  in  an  axially  inhomogeneous  cylindrical  plasma-vacuum  system.  .A  new  method 
of  solving  this  complicated  problem  is  proposed  and  carried  out. 

.An  analysis  of  ICRF  propagation  in  an  axially-inhomogeneous  cylindrical  plasma-vacuum 
system  is  undertaken.  Due  to  the  conservation  of  the  flux  embedded  in  the  plasma,  the  radius 
of  the  plasma  cylinder  changes  with  the  v-arying  axial  magnetic  field  as 

i  1C  RH2)B[z)  =  fijBo  =  R{z)  = 

Because  of  this  nonconstant  plasma- vacuum  interface  R{z),  matching  boundary  condi¬ 
tions  across  the  plasma-vacuum  interface  poses  a  major  problem.  The  traditional  (r,  z) 
cylindrical  geometry  is  transformed  into  a  new  geometry  (p,0  which  is  selected  to  facili¬ 
tate  the  matching  of  the  boundary  condition  and  to  simplify  the  form  of  the  plasma  field 
equations.  The  new  vzuiables  p  and  4  are  defined  as 

’’  -  R(z)  ’  -  R(z) 

In  this  new  coordinate  system,  the  boundaries  are  located  at 

•  p  =  1  :  Plasma- vacuum  interface. 

•  p  ~  K  :  Conducting  wall. 

This  new  coordinate  system  (p,^)  has  a  much  simpler  boundary  at  the  plasma- vacuum 
boundary  (p  =  1)  as  opposed  to  the  varying  radius  of  the  plasma- vacuum  interface  r  —  R{z) 
of  the  (r,  2)  coordinate  system.  The  new  (p,0  geometry  includes  a  conducting  wall  with  a 
radius  that  changes  in  proportion  to  the  plasma  column  radius.  This  is  not  an  unrealistic 
assumption  since  the  v^acuum  chamber  necks  down  near  the  B-field  maximum  in  most  mirror 
machines  in  order  to  locate  the  magnets  closer  to  the  plasma. 
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Actual  Experimentai  Geometry  of  Mirror  Machine 


Figure  4:  Geometry  of  the  Cylindrical  Plasma- \  acuum  System.  An  experimental  geometry 
(The  MIT  Tandem  Mirror  Rocket)  is  shown  schematically  at  top,  the  model  geometry  used 
in  the  analytical  theory  is  shown  at  bottom.  The  source  antenna  is  a  dual  half-turn  loop  for 
the  MIT  PPEX  experiment  and  a  dipole  loop  for  the  model. 
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The  steps  taken  in  the  solution  of  the  KM  fields  in  this  inhomogeneous  medium  are  as 
follows: 

1.  Solve  for  simrce  fields  generated  by  dipole  loop  in  vacuunn  A  complete  solution,  i.e. 
both  coefficient  and  functional  dependence  of  the  fields  is  possible  because  the  source 
term  is  known. 

2.  Solve  ff)r  the  full  vacuum  fields 

3.  Starting  from  Maxwell’s  equations  and  assume  the  following 

(a)  m  =  0.  Only  axisymmetric  mode  is  considered. 

(b)  Etf,  Er-  This  assumes  that  the  field  is  of  order  e  lower  than  Eg  or  Er, 
where  e  =  mgjmi. 

(c)  Constant  radial  density  profile  (skin  current  model).  .A  macroscopic  (non-RF) 
surface  current  is  assumed  to  exist  on  the  plasma- xTicuum  boundzu'y.  This  sur¬ 
face  current  generates  an  MHD  pressure  balance  [14]  that  can  support  a  radially 
c<»nstant  density. 

(d)  Slow-z  approximation.  The  magnetic  field  is  assumed  to  be  a  slowly  varying 

function  of  z.  This  means  that  fcy  --  1/a  »  I /Lb  .  where  ky  is  the  transverse 

wave  number  of  the  ICRF  waves,  fc*  is  the  axial  wave  number,  a  the  radius  of  the 
plasma,  and  and  Lb  is  the  axial  scale  length  of  the  ms^netic  field.  This  assumption 
simplifies  the  cross-coupling  terms  in  p  and  ^  in  the  differential  equations. 

From  these  assumptions,  one  can  derive  two  sets  of  coupled  differential  equations  link¬ 
ing  d>p  and  Xp,  where  Eg  =  iii/<j>p/r,  Er  =  iuJXp/r,  and  -  iu)<i>'p/r  are  the  RF  electric 
and  magnetic  fields  in  plasma. 

+ kimmoxp + =  0  (5) 

= 0  (6) 

where  R{^)  is  the  plasma  radius  as  a  function  of  ^  and  Ki(0  and  £2(0  elements 
of  the  dielectric  tensor  K  where 


VxE  =  ~—KE  ; 
c*  — 


K  = 


K,  K2  0 
-K2  A'l  0 
0  0  £3 
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1.  The  roeffii  ients  of  the  differential  equation  K\  2(0^*(0  decomposed  into  three 

parts;  its  value  as  ^  •  x.  its  resojiance  pole  contribution,  and  a  smooth  |)ortion,  i.e. 


KM)R^{i)-S^  ■  5.^(0 


5,^(0 


resonance  pole  smooth  portion 


Likewise 


KM)R\i.)  =  -i  +  D^iO 


resonance  pole  smooth  portion 


5.  The  functional  dependence  of  the  Resonance  Pole  term  is  selected  to  be 


Ti^tub  — 


where 


(1  +  ecn)/o^  +  «(1  +  cijc)  (j  +  ecuifU  +  +  ^22^) 


/o  =  m  =  Ha) 


.«) = (i 


where  Ra  is  the  radius  of  the  resonance  layer  of  plasma.  Note  that  has  two 
components.  The  first  component  is  for  the  subtraction  of  the  pole  at  ^  while 

the  second  component  treats  the  pole  at  ^  =  — $q. 

6.  Assume  a  Bessel-type  radial  solution  and  perform  a  Fourier  transform  on  <i>. 

2. 

Solving  this  yields  the  functional  form  for 

<t>{pjk^)  =  <f>{a,k^)pJi{ap) 

7.  Boundary  condition.  Application  of  the  continuity  of  fields  at  the  plasma- vacuum 
boundary  (p  =  1)  yields  the  following  for  0(Q,fc{). 


i  /  I.  x  1  [  5  ^ 

aJoia)  -  \yP  ^  . 
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where 


F{p,k^) 

6’(p,fc«) 


"  (k)  '  ^u„re(P,fc€) 


(8) 

(9) 

(10) 


It  is  seen  that  (l>  is  completely  determined  from  Equation  7  since  the  source  fields  ^uurr 
is  completely  known. 

8.  Perform  a  Fourier  transform  in  (  on  the  two  coupled  differential  equations  in  4>p  and 
Xp.  The  Fourier  transform  of  the  resonance  pole  5,^^  and  D^jb  are  performed  using 
residue  calculus  to  transform  the  analytically.  The  smooth  term  Sm  and  Dtrn  are 
transformed  first  by  curve  fitting  them  with  a  given  function  in  minimizing  the  error 
in  the  curve  fit  by  adjusting  some  coefficients,  then  Fourier  transforming  the  curve 
fitting  function  analytic2dly. 

The  resulting  equations  contain  convolution  integrals  because  of  the  ^  dependence  of 
the  coefficients  Ssm-,  and  Z>^. 


’ 

klXp{a,k^)  ^  feo  I 


^  f  [S^kf  -  k^)  r  S,^(k(  -  A')]  i^(a,  k^)dk' 
-if[£f.m(k(  -  k^)  +  £>sub(k(  -  ^p(a,k')dk' 

-S^Xpia.k() 


}  =  0 


[  -iDoo<t>p{a,ki) 


{-kl  -  a^)^j,{a,k()  +  14  I 


[  +  /  -  k'^)  +  S,^{ki  -  i^{a,  k'^)dk'^ 

+*/  -  k')  +  D^(k(  -  fc')]  ip(a,k')dk'  I  ^  ^ 

+  S*,0p(Q,fe«) 

[  +iDocXp{a,h) 


9.  Finite  differencing  the  two  coupled  differential  equations  yields 


khT  +  ^ 


^Sj^ixT^ki-i'£Dj.i<i>r^k  +t^{SooXT-iD^<t>^)=0  (11) 


(fcj  f  a,^)<^7  f  fco^ 


L  (  I 


+  (s»«7  +  WocxT)  =  0 

(12) 
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where 

kj  -  k^  ;  k,  T  fc' 

xT  - 

<p^  =  4>^(a"'.k^) 

Si_i  =  [5,„(Ar(  -  k[)  ^  S,^ik^  -  A:;)] 

-  k'^)  -  D,^ik^  -  Af^)] 

The  two  coupled  difference  equations  could  be  written  as  two  sets  of  matrix  equations. 

4(^)*X  =  =Si(a) 
il2(")  ■  X  ="  -^2(0) 

They  are  then  solved  simultaneously  to  yield  an  a  for  each  as  well  as  x-  The 
(*>  is  solved  directly  from  the  boundary  conditions.  These  are  then  inverse  Fourier 
transformed  to  return  the  ICRF  fields  £„  and  B,. 

10.  A  figure  of  merit  is  selected  as  the  antenna  impedence.  Rpi.  This  is  found  by 

E  J‘dV 

which  is  the  volume  integrated  absorbed  power.  This  figure  of  merit  can  then  be 
maximized  for  different  field  and  plaisma  conditions. 

■A  code,  ('VLWAVE,  is  written  in  order  to  implement  the  theoretical  analysis  detailed 
above.  The  computed  ICRF  fi,  from  CYLWAVE  for  a  resonance-in-chamber  case  is  plotted 
on  Figure  5.  It  is  seen  that  the  ICRF  waves  decay  very  quickly  after  leaving  the  launching 
zmtenna.  For  a  Resonance-out-of-chamber  case  seen  on  Figure  6,  it  is  seen  that  the  absence 
of  a  resonance  results  in  more  oscillations  and  larger  field  amplitudes  after  the  ICRF  waves 
are  launched  from  the  anteima. 

A  mirror  ratio  scan  is  performed  in  order  to  calculate  the  mirror  ratio  for  the  optimum 
wave-plasma  coupling.  The  plasma  impedence  Rpi  at  different  mirror  ratios  is  plotted  on 
Figure  7.  The  behavior  of  the  Rpi  as  a  function  of  the  mirror  ratio  is  in  good  general 
agreement  with  the  resonance  volume  theory,  which  simply  states  that  the  dissipation  is 
proportional  to  the  volume  of  plasma  that  is  close  to  resonance. 
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The  behavior  of  the  antenna-plasma  coupling  inipedence  as  a  function  of  resonance  lo- 
t  at  ion  is  also  perf<»rmed.  The  rt'siilting  Rpt  versus  is  plotted  on  Figure  8.  The  coupling 
iinpedence  is  seen  to  peak  in  two  locations  -  rescmance  at  plasma  center,  and  resonance 
directly  iinderneath  the  antenna.  In  fact,  the  resonance-under-antenna  case  shows  peak  in 
plasma  loading  impedence  that  is  just  as  large  as  the  resonance- at -center  case.  This  indicates 
that,  depending  on  the  mirror  ratio,  locating  the  resonance  close  to  the  antenna  may  be  just 
as  important  for  coupling  as  locating  the  resonance  at  the  plasma  center.  Thus  an  optimum 
coupling  condition,  from  a  wave- propagation  point  of  view",  would  be  to  have  the  resonance 
located  either  directly  underneath  the  antenna  or  at  the  center  of  the  machine. 

It  must  be  noted  here  that  the  analytical  study  described  here  hais  two  significant  adran- 
tages  over  other  studies.  First,  the  axial  inhomogeneity  and  cylindrical  geometry  are  both 
treated  in  a  consistent  manner  in  the  analysis,  and  secondly,  the  subtraction  of  the  resonance 
pole  performed  in  this  study  allows  the  ICRF  waves  to  propagate  up  to  and  through  the 
resonance.  These  advantages  demonstrate  the  power  of  this  method  of  analysis. 

.\s  summary  of  this  analysis  w'e  can  draw  the  following  conclusion:  A  high  coupling 
efficiency  can  be  obtained  if:  (1)  The  magnetic  filed  is  made  as  uniform  as  possible;  (2)  The 
rf  frequency  is  in  resonance  with  ion  cyclotron  frequency  at  minimum  field;  and  (3)  The 
coupling  efficiency  increases  with  plasma  radius.  This  w'ill  lead  our  future  experiment. 
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Hgure  o:  K’KF  pl<»l  as  function  of  Ip.i).  Waves  are  launched  from  dipole  ItJop  antenna 
located  at  p  -  1.397,  ^  -  o.Od.  For  this  case.  K’RF  resonance  location  is  h)cated  l.\ 

the  chamber  at  ^  -  i().977.  65  spectral  points  were  taken  to  generate  this  prohle.  The 
geometrv  is  /?(,  0.77(f.  -  n..5u.  where  a  is  the  radius  of  the  wire  loon  antenna. 
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Figure  7:  Plasma  Impedence  Rpi  (Solid  line)  and  Normalized  Resonance- volume  (Open  Circle 
Line)  calculations  performed  for  different  mirror  ratios.  Plasma  parameters  are:  n,  =  1.0  x 
10'‘cm  ^re  -  lOeV. 
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Figure  8:  Plasma  Impedence  Rpi  versus  rest)nance  location  ^res-  Plasma  parameters  are: 
n,  =  1,0  X  10^‘cm'^,  T,  =  lOeV.  Note  that  the  peaks  in  Rpi  are  located  at  ^  =  0  field 
minimum,  and  ^  ■  5.00  antenna  location. 
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3.2  Computational  Study 

Although  the  analytical  treatment  «>f  the  previous  section  yielded  intuitive  results  of  the 
dependence  of  on  plasma  parameters,  it  nevertheless  is  too  simplified  a  treatment  to 
apply  to  a  real  system.  I'he  main  simplifications  that  were  made  in  the  analytical  treatment 
was  the  use  of  only  one  eigenmode  and  homogeneous  magnetic  field.  A  typical  solution  that 
proved  to  be  useful  requires  the  solution  of  about  30-50  modes  which,  upon  summation, 
form  the  resjilting  field  solutions. 

The  goal  of  this  computational  analysis  is  to  provide  a  better  modeling  of  the  wave- 
plasma  coupling  problem  than  the  simplified  results  of  the  analytical  solution.  However,  this 
computational  analysis  is  not  simply  a  purely  numerical  simulation.  Much  analytical  work 
is  done  to  reduce  the  fields  to  functional  form  ,  i.e.  In{krr),  K„{krr),  J„{kxir)  before  the 
numerical  “number-crunching”  is  handed  over  to  the  codes. 

The  BEACH  code  extends  the  analytical  approach  taken  in  the  cylindrical  geometry  so¬ 
lution  of  the  previous  section.  It  solves  for  the  ICRF  wave  propagation  in  an  inhomogeneous 
axial  magnetic  field.  The  inhomogeneous  axial  geometry  is  handled  by  separating  the  axial 
length  into  discrete  axial  “slices”,  each  of  which  contains  a  homogeneous  axial  magnetic  field. 

The  code  first  finds  the  eigenmodes  that  exist  in  each  axial  slice  by  satisfying  both  the 
infinite  plasma  dispersion  relation  and  the  boundary  conditions.  The  boundary  condition  of 
the  plasma-  vacuum  system  can  be  written  as; 


ciHp  ^  C2Ep  =  eiJ^  (13) 

csHp  *  c^Ep  =  e2J^  (14) 

The  solution  f<»  the  above  are; 

Hp  =  (15) 

Ep  =  (16) 

D~CiC4-C2C3  (17) 

The  term  D  is  then  the  boundary  condition  of  the  bounded  plasma-vacuum  system.  Solving 
this  boundary  condition  coupled  with  the  dispersion  relation 

Ski  kl[kl{S  +  P)-  k^{SP  -t-  RL)]  +  P{kl  -  klR){kl  -  k^L)  =  0  (18) 

yields  the  eigenmodes  of  the  cylindrical  system. 

After  solving  for  the  eigenmodes  of  the  system,  the  code  then  computes  the  excitation 
coefficients  of  the  initial  axial  slice  as  driven  by  the  source  term  (Antenna).  The  excitation 
coefficient  of  the  next  axial  slice  can  then  be  determined  by  imposing  two  conditions  of 
continuity  on  the  boundary  betw’een  the  two  axial  slices;  conservation  of  the  integrated 
axially  propagating  Poynting  Flux  f  S2rdr  ,  and  the  continuity  of  wave  phase.  These  two 
conditions  uniquely  determine  the  excitation  coefficients  of  the  different  eigenmodes  in  the 
next  axial  slice.  This  analysis  is  then  performed  for  all  axial  slices. 

Computational  analysis  was  performed  on  two  types  of  plasma  regimes;collisionless  cold 
Plasma,  and  collisional  cold  plasma.  Different  physics  are  inherent  in  each  different  regime. 
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3.2.1  CoUisionless 

Initially,  the  calculations  of  the  IC'HF  field  propagation  were  performed  using  a  cold,  colli¬ 
sionless  plasma  model.  This  ensures  that  there  exist  no  complex  roots  that  could  damp 
out  the  axially  propagating  Poynting  Flux,  thus  enabling  a  good  check  of  the  model's  ability 
to  conserve  power  and  basic  phsical  insight.  Because  there  is  no  dissipation  mechanism  (Cold 
plasma,  no  collisions),  some  effects  are  expected  physically  from  this  simple  propagation. 

.Several  different  antenna  geometries  are  used  to  launch  different  modes.  The  two  antenna 
geometries  that  were  studied  were  the  FULL  TURN  LOOP(FTL)  and  DUAL  H.ALF  TUR.N 
LOOP(DHTL).  The  baseline  model  is  driven  by  a  Full  Turn  Loop  (FTL)  antenna.  Because 
it  is  symmetric  in  6,  the  full  turn  loop  antenna  excites  only  the  n=0  mode. 

The  solution  of  the  dual  half  turn  loop  (DHTL)  geometry  poses  a  more  difficult  problem 
than  the  full  turn  loop  (FTL)  solution.  The  main  difference  betwwn  the  DHTL  and  the 
FTL  lies  in  the  J(fcj)  driving  term,  specifically  the  radial  feeders.  This  is  different  from  the 
Full  Turn  Loop  (FTL)  geometry  in  which  the  contribution  from  the  radial  feeders  cancel 
out  (Since  one  leg  is  entering  and  the  other  is  exiting,  net  is  zero).  The  modelling  of  the 
driving  term  J(A:j)  for  the  DHTL  is  detailed  in  a  full  report  in  preparation. 

The  wave  propagation  of  both  type  of  antenna  is  quite  similar.  Figures  9  and  10  show 
the  propagation  characteristics  of  the  ICRF  and  fields  eis  the  waves  propagate  from 
the  launch  point  (z=0)  towards  the  resonance  at  mid-plane.  The  results  that  are  readily 
obser^Tible  from  this  coUisionlcss  propagation  of  the  n-0  mode  towards  resonance  arc  given 
below: 


1 ,  IB;  Increases  -  This  is  an  expected  result.  Stix  shows  (Stix,  Theory  of  Plasma  Waves, 
p.5l.  Eq(18)  )_15  that  for  a  homogeneous  lossless  medium,  the  group  velocity  is 


V 


9  ~ 


IVo 


(19) 


When  resonance  is  approached,  the  group  velocity  slows  down,  resulting  in  an  increased 
energj'  density  Wq.  This  increased  energy  density  is  apparent  in  the  increase  in  \B.\, 
since  the  energy  density  \'aries  as 


Ho  = 


+  to 


(20) 


2.  Increase  in  axial  and  radial  k*  -  As  resonance  (u;  =  ui^)  is  appraoched,  the  left 
hand  cut-off  (fe|  =  fc^Z,)  increases,  causing  am  increase  in  the  resulting  k;  eigenmodes. 
This  is  consistent  with  the  slowing  down  of  the  group  and  phase  velocities  as  resonance 
is  approached. 

3.  High  edge  electric  fields  £+ ,  £_  and  E,  that  propagate  inward  -  This  is  a  result 
that  could  be  consistent  with  experimental  obser\'ation.  Golovato  points  out  (Personal 
communications)  that  although  a  high  edge  electric  field  is  predicted  by  the  ANTEN.4 
code,  experimentally  no  edge  electric  field  is  detectable.  The  BEACH  simulation  shows 
that  the  reason  that  no  edge  electric  field  is  detectable  is  that  the  edge  spike  moves 
radially  inward  into  the  plasma  as  resonance  is  approached. 
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N*0.  NB*  1.50«10“.  W/VCIO*  0.788  .  COLD  COLUSIONLESS  MODEL 
NO.  MODES  >  00 


Figure  9:  ICRF  Propagation  in  a  Cold  ColHsionless  Plasma,  Full  Turn  Loop  Antenna, 
u;/u;ci0  =  0.788.  Waves  are  launched  by  the  ICRF  antenna  at  2  =  0cm  and  propagate 
towards  resonance  at  2  —  50cm. 
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N«0,  NE«  1  60*10**  W/VCIO*  0.788  .  COLD  COLLISIONLESS  MODEL 
NO.  MODES  -  80 


Figure  10:  ICRF  Ez  Propagation  in  a  Cold  Collisionless  Plasma,  Full  Turn  Loop 
ijj/u!ciO  —  0.788 


Antenna, 
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1.  Higher  order  modes  excited  at  higher  amplitudes  N  for  lower  initial  - 
This  result  reflects  the  effect  of  the  antenna  excitation  profile  (Jeii^z))  *^*1  excited 
EM  field  coefficients.  For  a  lower  the  value  of  the  left-handed  cut-off  in  the 

infinite  plasma  dispersion  relation  k^L  decreases.  Since  the  Jeik.)  profile  remains 
constant,  this  means  that  .ALL  modes  are  now  excited  to  higher  values. 


3.2.2  Collisional 

The  presence  of  collisions,  of  course,  implies  a  finite  temperature.  The  Cold  Plasma  Colli¬ 
sional  Model  means  that  the  expansion  of  the  plasma  Z-function  excludes  temperature  effects 
and  keeps  only  terms  of  0(i/,7u;).  The  code  results  for  the  Ceise  of  ICRF  propagation  in  a 
cold  collisional  plasma  is  shown  in  Figure  11. 

The  presence  of  collisions  manifests  itself  in  a  complex  dielectric  tensor. 


where 


c 


5  -iD  0 
iD  S  0 
OOP 


(21) 


5=i{P  +  I)  (22) 

D=\{R-L)  (23) 

The  presence  of  collisions  in  a  cold  plasma  produces  the  following  expansions  for  S,  D,  and 
P 


w, 


pa 


-  u;^  +  i/f) 


^f) 


(24) 


a  u»(u>-)-tl'ca) 


(25) 

(26) 


Using  these  complex  elements  of  the  dielectric  tensor,  the  dispersion  relation  can  be  solved 
with  the  geometric  boundary  condition  to  yield  a  set  of  eigenmode  solutions  for  this  plasma- 
x'acuum  geometrj’.  The  difference  between  these  solutions  and  those  of  the  cold  collisionless 
plasma  is  that  the  eigenmodes  are  now  complex  -  these  waves  are  now  damped  as  they 
propagate  towards  the  cyclotron  resonance.  This  damping  transfers  energv-  from  the  waves 
to  the  plasma,  causing  plasma  heating. 
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Figure  12:  ICRF 
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4  Thrust  Balance  Development 

I'hfre  is  a  need  Id  develop  an  accurate  and  sensitive  target  balance  to  measure  the  momen¬ 
tum  transfer  in  a  pla>ma  system,  such  as  an  electric  propulsion  device,  plasma  propulsion 
systems  and  in  the  divertor  of  a  tokamak  fusion  reactor.  In  electric  propulsion  research  the 
common  method  to  measure  the  thrust  is  a  balance  consisting  of  the  suspended  thruster  and 
a  displacement  sensor.  The  thrust  is  determined  from  the  displacement  of  the-thruster  when 
it  recoils  at  the  firing  of  the  thruster.  This  method  has  several  drawbacks:  The  thruster  is 
heavy.  It  is  tied  to  heavy  electric  power  cables  and  gas  feed  line  and  therefore  has  a  very  high 
motion  resistance.  There  are  frictional  losses  due  to  the  suspension  of  the  heavy  thruster 
and  the  transfer  of  displacement  through  pulleys  to  the  sensor.  In  the  future,  the  header 
will  probably  need  to  be  cooled  for  high  power  operation.  The  motion  resistance  will  be 
much  higher  with  added  cooling  lines  and  heavier  power  cables.  These  factors  will  reduce 
the  sensitivity  and  accuracy. 

For  the  purpose  of  studying  the  momentum  transfer  in  a  fusion  plasma  system  there 
are  two  methods  being  used  by  other  groups.  One  is  to  attach  the  target  through  a  long 
shaft  to  the  diaphragm  of  a  commercial  pressure  transducer.  The  shaft  is  suspended  by 
two  strings.  The  suspended  target  and  the  shaft  system  behaves  like  a  pendulum  and  the 
pressure  transducer  serves  as  a  displacement  sensor.  The  other  method  is  similar  to  the 
suspended  target  and  shaft  assembly  whereas  the  sensor  consists  of  a  magnet  and  a  pick-up 
coil.  The  force  can  be  determined  from  the  voltage  induced  on  the  pick-up  coil  from  the 
mettion  of  the  magnet.  The  sensor  unit  has  to  be  housed  in  a  magnetically  shielded  box. 
The  l«»ng  arm  of  these  two  systems  will  make  them  difficult  to  be  installed  in  the  plasma 
stream  in  a  closed  system  like  tokamak.  None  of  the  systems  mentioned  above  has  provided 
adequate  means  to  deal  with  the  problem  of  mecheuiical  vibrations,  electric  and  magnetic 
interferences. 

The  target  balance  system  developed  in  this  work  consists  of  two  identical  plates  sus¬ 
pended  by  needle  points.  The  device  is  compact,  light  and  is  symmetric  with  respect  to  the 
equilibrium  position  as  either  or  both  plates  is  impinged  by  plasma  from  opposite  directions. 
The  system  can,  therefore,  be  inserted  in  the  plasma  stream  in  a  closed  system  like  tokamak. 
The  mechanical  vibrations  can  be  eliminated  as  a  common  mode  to  the  two  identical  moving 
pendulum.  The  electric  and  magnetic  pick-ups  have  been  reduced  to  a  negligible  level  by  the 
use  of  differential  amplifiers,  shieldings  and  dynamic  bandpass  filters.  The  detailed  devel- 
opuiental  process  of  the  pendulum  system,  electronic  circuit,  the  analysis  and  measurement 
are  described  in  the  following  subsections. 
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4.1  The  Balance  Development 

The  purpose  of  this  work  is  to  find  a  suitable  method  to  measure  the  thrust  in  the  Tandem 
Mirror  Plasma  propulsion  device  Ij  2i  ,'3;  [16  .  At  present  the  device  is  operated  in  the  low- 
power  level  and  at  a  high  specific  impulse  (  about  I3,t)00  sec  ).  The  thrust  le\el  is  very  low 
and  is  estimated  to  be  at  about  100  mN.  To  measure  such  a  low  thrust  a  sensor  must  have 
the  sensitivity  and  accuracy  of  better  than  .1  mN.  The  balance  has  to  be  very  compact  so 
that  it  can  be  mounted  in  the  limited  space  of  the  exhaust  chamber. 

The  pr«»pulsion  system  is  presently  operated  in  a  pulsed  mode  and  plasma  is  created  by 
microwave  radiation  and  heated  to  high  temperature  by  rf  power.  The  magnetic  impulse 
and  mechanical  vibrations,  rf  pick-ups,  electric  field,  electromagnetic  noises  from  the  power 
supply  systems  are  many  orders  of  magnitude  stronger  than  the  anticipated  thrust  level. 
Undoubtedly  they  must  be  eliminated.  At  the  beginning  we  were  hoping  to  adapt  the  known 
methods.  As  described  above  none  of  them  are  adequate  for  our  system.  The  possibility  of 
using  torsion  wire  and  strain  gauge  has  also  been  carefully  examined  and  found  not  desirable 
either. 

In  order  to  keep  the  system  as  simple  and  compact  as  possible,  we  have  given  the  simple 
pendulum  a  serious  consideration.  However,  it  loses  its  simplicity  and  compactness  when 
the  position  sensing  system  is  incorporated.  After  weighing  all  the  odds  we  decided  a  system 
consisting  of  double  plates  is  our  best  bet.  We  were  encouraged  by  the  initial  success  of 
a  pendulum  and  fixed  plate  system  as  shown  in  Figure  13a.  Both  the  moving  and  fixed 
plates  are  identical  in  size  and  shape  and  are  made  of  aluminium.  The  displacement  of 
the  pendulum  is  translated  into  variation  of  capacitance  of  the  plates.  The  capacitance 
rneastiring  method  will  be  described  later.  The  signals  for  the  balance  in  the  air  and  in  the 
vacuum  chamber  are  shown  in  Figures  14a  aiid  14b  respectively.  The  oscillations  in  the  air 
were  produced  by  mechanical  shock.  The  oscillation  in  the  vacuum  is  produced  by  pulsing 
the  magnetic  field.  As  can  be  seen  from  these  figures,  the  signal  taken  in  the  air  damped 
quickly  whereas  there  was  almost  no  damping  in  the  \'acuuin.  This  shows  that  the  friction 
loss  at  the  needle  support  is  negligible.  The  magnitude  of  the  magnetic  impulse  is  very  large 
and  can  be  eliminated  by  dividing  the  plates  into  very  thin  strips. 

The  mechanical  vibrations  are  the  second  order  noise  which  appears  when  the  magnetic 
impulse  is  removed.  As  shown  in  Figure  14c,  such  noise  is  random  emd  could  not  be  reduced 
with  suspension  methods.  To  eliminate  the  mechanical  noise,  both  plates  are  suspended  and 
swing  as  an  independent  pendulum  (  Figure  13b  ).  Since  both  plates  respond  identically 
to  the  noise,  no  net  change  in  the  capacitance  is  produced  and  therefore  no  net  signal  is 
recorded  from  the  noise. 

In  the  plasma  stream,  the  conducting  plate  acts  as  a  big  Langmuir  probe  which  will 
detect  the  electric  field  and  current.  These  electric  field  are  the  third  order  effects  which 
appears  when  the  vibrational  noise  is  cancelled  out.  As  shown  in  Figure  14d,  these  electric 
fields  are  still  very  large.  To  reduce  the  electric  field,  insulating  materials  such  as  ceramic  or 
silicon  are  used  on  the  surface  facing  the  plasma.  Further,  a  set  of  plates  with  conducting 
strips  are  placed  above  the  support.  The  electromagnetic  noise  from  the  inverse  pendulums 
can  then  be  used  to  cancel  the  noise  picked  up  by  the  regular  pendulums  below  the  support. 
The  inverse  pendulum  above  the  support  swings  opposite  to  the  regular  pendulum  below  so 
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balance  consisting  of  a  simple  pendiil 


igiire  J4:  Signals  from  the  thrust  balance:  (a)  A  damped  oscillation  from  a  mechanical 
lock  in  the  air,  (b)  a  nundamping  oscillation  produced  by  magnetic  impulse  in  the  vacuum, 
)  mechanical  noises  and  (d)  electric  must's. 


that  the  signs  of  the  voltages  from  these  two  sets  «>f  capacitor  plates  are  opposite,  whereas 
the  ><igns  of  pirk-ii|)  voltages  from  electric  and  magnetic  noise  are  the  -ame.  Ry  using  a 
rlifferential  amplifier  the  pick-up  signals  are  cancelled  out  and  the  real  signals  are  summed. 
This  system  is  shown  in  Figure  In. 

Since  the  inverse  pendulums  arc  smaller  than  the  regular  ones,  the  noise  can  not  be 
removed  completely.  To  further  avoid  the  electric  field  effect  due  to  plasma,  the  plate  in 
contact  with  the  plasma  is  made  entirely  with  an  insulator,  i.e.,  no  conductor  strips  even  on 
the  back  face.  As  depicted  in  Figure  16,  the  device  is  now  a  quadruple  pendulum  system 
where  there  is  a  set  of  double  capacitor  plates  above  the  support  and  another  set  of  double 
capacitor  plates  below  the  impact  plates  which  are  now  made  of  insulators.  The  double  ca¬ 
pacitor  plates  consist  of  conductive  metal  strips  used  to  monitor  the  signals.  These  capacit<jr 
plates  are  enclosed  in  a  conducting  box  to  shield  out  the  electric  and  magnetic  pick-up.  .Any 
residue  pick-up  can  be  further  eliminated  with  the  use  of  differential  amplifiers.  This  method 
effectively  reduces  the  pick-up  and  enhances  the  sensitivity.  The  magnetic  noise  produced 
by  the  power  supply  was  the  most  difficult  to  deal  with  and  can  not  be  eliminated  by  the 
means  described  above.  It  was  reduced  with  a  carefully  designed  active  filter. 

4.2  The  Signal  Treatment 

The  simplified  electronic  circuit  diagram  is  shown  in  Figure  17.  To  detect  the  small  change 
in  capacitance  of  the  two  pendulums,  a  carrier  radio  frequency  V(^ltage  at  30  kHz  is  applied 
to  plates.  There  are  two  pairs  of  signals,  one  pair  from  the  top  capacitor  plates  and  the  other 
from  the  bottom  capacitor  plates.  The  signals  from  the  each  pair  of  plates  are  carried  by 
a  pair  of  twisted  wires  to  the  inputs  of  a  differential  amplifier.  Since  the  inverse  pendulum 
plates  on  the  top  swing  in  the  opposite  direction  to  the  pendulum  plates  at  the  bottom,  The 
signs  of  the  signal  from  the  output  of  the  differential  amplifiers  are  opposite.  These  outputs 
are  again  summed  by  the  third  differential  amplifier.  Therefore  the  overall  signal  is  enhanced 
after  all  the  noises  are  eliminated. 

The  carrier  frequency  is  removed  from  the  signal  by  rectifiers.  The  residual  noise  is 
practically  cancelled  by  the  differential  amplifiers.  The  magnetic  ripples  from  the  firing  of 
SCRs  in  the  power  supplies  are  the  most  stubborn  problem.  The  frequency  of  this  noise  has 
fallen  in  the  range  of  the  carrier  frequency.  This  noise  is  filtered  by  a  very  carefully  designed 
active  b2mdpass  filter.  A  high  quality  clean  signal  can  finally  be  obtained. 
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Figure  J6:  The  working  version  of  the  thrust  balance  consisting  of  noncondiicliiig  largel 
plates  and  sensor  capacitor  plates  with  fine  conducting  strips  encKised  in  shielded  Imxes  on 
the  lop  and  Imttom. 
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NOTE: 

OS:  oscillator  40KHz,100V 
P1,P2:  double  pendulum 
A1,A2:  pre-amplifier 
F1.F2:  active  filter 
REi.RE2:  rectifier  it  filter 
A3:  differential  amplifier 
A4:  output  amplifier 
F3:  anti-DC-dri f t  filter 
F4:  low  pass  filter 


Figure  17:  Block  diagram  of  the  electronic  circuit 
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4.3  Experimental 

The  two  pendulums  of  the  thrust  balance  are  not  exactly  identical  and,  therefore,  their 
natural  frequencies  are  not  equal.  When  a  mechanical  impulse  is  applied,  a  beat  as  shown 
in  Figure  18a  is  produced.  The  frequency  of  the  beat  is  about  one  twentieth  of  the  natural 
frequency  or  the  period  is  20  times  longer.  This  beat  frequency  will  become  smaller  when 
two  natural  frequencies  are  made  to  be  closer.  A  great  effort  has  been  made  to  tune  the 
two  pendulums  to  make  their  natural  frequencies  match.  H.nvever,  a  perfect  match  can  not 
be  achieved  and  it  can  only  make  the  beat  period  longer.  The  beat  can  not  be  eliminated 
completely.  Currently  the  thruster  is  operated  in  a  pulsed  mode  and  is  fired  in  the  first  half 
cycle  of  the  natural  oscillation.  The  trace  of  the  first  cycle  in  Figure  18a  is  expended  zis  is 
shown  in  Figure  18c.  This  expended  trace  shows  that  the  noise  level  is  negligible.  Figure  18d 
presents  the  relationship  between  the  plasma  density  and  temperature  traces  from  a  triple 
probe  and  the  pendulum  signal.  This  shot  was  made  using  the  pendulum  system  shown  in 
Figure  15  where  the  signal  died  out  after  one  cycle.  It  should  continue  to  oscillate  without 
damping  and  therefore  it  may  not  be  the  valid  signal.  In  order  to  determine  the  validity  of 
the  measurement  it  is  necessary  to  carry  out  detailed  analysis  which  w’ill  be  given  in  the  next 
section.  An  exhaustive  study  of  the  residual  noise  from  the  power  supply  reveals  that  they 
can  sometime  cause  misleading  false  signal  even  at  a  very  low  level.  The  complete  elimination 
of  such  noise  is  necessary  and  was  accomplished  with  the  use  of  dynamic  bandpass  filters. 
Finally  Figure  19  presents  the  signal  obtained  with  the  working  pendulum  system  shown 
in  Figure  16  using  the  electronic  circuit  shown  in  Figure  17.  The  signed  continued  to  oscillate 
without  damping  but  W'as  not  symmetric  with  respect  to  the  baseline.  This  requires  a  solid 
explanation  in  order  to  convince  cuirself  that  it  is  in  fact  a  correct  signal.  The  amplitude  of 
the  signal  rises  from  a  negative  value  to  a  constant  symmetric  oscillation  after  about  three 
cycles.  This  asymmetric  property  is  found  to  be  due  to  the  electric  circuit  as  demonstrated 
by  the  circuit  analysis  given  in  the  next  section. 


4.4  Analytical  Solutions 


As  discussed  above  in  order  to  V2didate  the  me2isurement  it  needs  to  carry  out  a  complete 
analysis  of  the  pendulum  system  and  electronic  circuit.  Let  us  wrrite  down  the  equations  of 
motion  for  a  simple  pendulum  as  showm  in  Figure  20.  Assuming  9  small,  the  equation  of 
motion  is 


<P9 

dt^ 


=  Ft  ~Ff-F, 


(27) 
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PENDULUM  CHARACTERISTICS 
ANALYSES 


Assuming  B  small,  the  equation  of  motion  is 


rrtl-^  =  Ft  -  Ff  -  Fg 


where  Ft  =  Thrust 
Fg  =  mgd 

A  a  damping  factor 
Let  T  ~  puke  length 

To=y/i/9 


Figure  20:  Simple  pendulum. 
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whore 


Ft  -  Thrust 

F'J  - 

Fg  =  mgB 

.4  =  adainping  fact«.)r 
Let  r  —  pulse  length 
To  = 


This  equation  has  been  solved  for  other  systems.  To  understand  the  characteristics  of  this 
system  the  procedure  for  solving  the  equation  of  motion  and  results  are  presented  in  the 
following. 

by  apply  Laplace  Transformation 


0(5)  = 


VFr, 


{S^ml  —  5.4  +  mg) 


(28) 


Examine  the  solutions  in  two  regimes: 

1.  r  <<  To 


fj'(t)  =  Fxr6(t),  a  8  function 

L.Ft  =  Fjt 


2.  r  >>  To 


Ft  =  Ft  ‘lit),  a  step  function 
L\Fr\  =  f 


In  each  regime  there  are  four  different  conditions: 

1.  /4  =  0,  no  damping 

2.  0  <  i4  <  Am2gl,  under  damping 

3.  A  =  4m2gl,  critical  damping 

4.  .4  >  4m2gl,  over  deunping 

In  T  «  To  regime,  apply  Inverse  Laplace  Transformation  to  the  equations  for  each 
condition,  we  obtain 


1.  For  /I  =  0 


^(0  = 


Ftt 

my/^ 


(29) 
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2  For  0  <  A  <"  \m2gl. 

m 

3.  For  .-4  =  4rr'*gl, 


‘2Ftt  .  -4  ,  .  ,  y/Aw?gl  .4^  . 

2m/'*""'  2m/-  ■'* 


^(<)  =  Vs/I') 


4.  For  .4  >  4m2gl, 


(30) 


(31) 


m 


These  solutions  are  plotted  in  Figure  21.  Comparing  the  analytical  solution  and  the  experi¬ 
mental  result  presented  in  Figure  19  it  appears  that  the  experimental  result  can  be  explained 
by  the  damped  solution.  But  it  has  shown  previously  that  the  friction  at  the  needle  support 
is  negligible.  The  damping  can  not  come  from  the  pressure  built-up  in  the  chamber  during 
the  shot  because  it  has  been  shown  that  there  is  no  observable  damping  in  the  chamber  at 
mtorr  pressure  range.  Therefore  this  is  not  a  true  result.  The  result  shown  in  Figure  19 
appears  to  be  real  and  in  agreement  with  the  undamped  solution  except  the  asynunetric 
property  in  the  first  three  cycle.  We  reason  that  this  is  an  effect  of  the  asymmetric  relative 
motion  of  the  capacitor  plates  and  the  electronic  circuit.  Because  the  relative  motion  of  the 
pair  of  plates  on  the  top  is  opposite  to  that  of  the  pair  «>f  the  bottom  plates,  the  asym¬ 
metric  effects  nearly  cancel  out  and  are  ignored  for  simplicity.  The  electronic  circuit  is  very 
complicated  and  there  are  high  nonlinear  elements.  Fortunately,  it  C2m  be  represented  by  a 
drastically  simplified  diagram  as  shown  in  Figure  22  and  two  simultaneous  linear  equations 
as  follows: 

^  ^  4.  Kl. 

~  Ti  ’ 

(32) 


( 


<u 


dt 


^  ^  J.  —  IOl 

dt  ^  T2  Ts’ 


where  Ti  =  R\C\  Juid  T2  =  R'2C2,  and  y  and  z  are  the  input  and  output  signal  respectively. 
The  solution  is  shown  in  Figure  23.  Except  for  a  slight  phase  different  the  analytical  solution 
closely  matches  the  experimental  data.  In  any  electric  system  phase  shift  is  common  and 
for  this  application  we  are  not  interested  in  its  value.  From  this  analysis  we  can  confidently 
claim  that  the  experimental  result  is  the  true  measurement  of  the  thrust  free  from  mechanical, 
electrical  and  magnetic  interferences. 
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PENDULUM  WITH  DIPPBRBNT  FRICTIONS 
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Figure  21:  I’he  waveforms  of  a  simple  pendulum  fc 


CIRCUIT  EQUATIONS 

THe  complete  drcuit  equation  is  high  order  nonlinear  and 
is  very  complicated.  However,  we  can  take  the  advantage  of 
the  short  time  constants  of  the  filters  FI  and  F2  and  that  the 
amplifiers  Al,  A2,  and  the  rectifiers  REl,  and  RE2  can  be 
considered  as  ideal  elements,  the  characteristics  of  the  circuit 
is  determined  mainly  fay  F3,  A4  and  F4.  The  circuit  can  then 
be  represented  by  the  diagram  and  two  simultaneous  linear 
equations  as  follows: 


t;-*,  c, 

Tj  «  Kj  Ct 


{ 


iSt  .  420.  10.. 

d*  “  dt  ^  iT* 
4jl  4.  ^  —  10. 

dt  ^  'T.  —  nr-  • 


Ta  -  Ta 

or  a  second  order  differential  equation  as 


_  <Pz  Tz.dz  z  dy 


Figure  22:  Simplified  circuit  diagram 


48 


CALCULATED  PENDULUM  SIGNAL 
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Figurt*  23:  Output  signal  wavpfurin  of  tlu-  simplifu-d  det  irii  (  im.i 


4.5  Experimental  Results  and  discussion 

I  he  initial  reMilts  obtained  using  this  balance  are  presented  in  Table  2.  The  plasma  tem¬ 
perature  and  density  were  monitored  at  a  radial  posititm  two  third  of  the  plasma  radius  by 
a  triple  probe.  The  temperature  at  the  center  is  monitored  by  a  spectrometer  using  Doppler 
broadening  method.  A  Langmuir  probe  biased  at  -100  is  used  to  measure  the  electron 
saturation  current. 

The  specific  impulse  is  12,852  sec  corresponding  to  the  plasma  temperature  of  172  e\  . 
The  measured  thrust  is  76  mN  in  agreement  with  prediction.  The  propulsive  efficiency  is 
68%.  This  can  be  considered  as  a  milestone  of  this  research. 

This  pendulum  is  made  in  the  laboratory  for  testing  purpose  and  is  first  of  its  kind. 
Therefore  there  is  room  for  improvement.  The  plates  are  made  from  printed  circuit  board 
which  consists  of  G-10  materials.  G-10  material  outgases  during  the  plasma  shot  which 
degrades  the  plasma  conditions.  The  new  plates  wdll  be  made  from  silicon  wafers.  For 
short  pulse  the  particles  can  all  be  absorbed  by  the  G-10  material.  For  long  pulse  the 
surface  reflection  has  to  be  kept  to  a  minimum  and  its  effect  has  to  be  factored  into  the 
measurement. 

There  are  two  major  advantages  of  this  balance:  (1)  Since  the  plates  are  not  parts  of 
propulsion  system  and  are  simple  and  inexpensive  to  build;  they  can  be  made  sacrificial  for 
high  power  operation.  (2)  The  balance  is  symmetric  in  both  directions.  It  is  particularly 
useful  to  study  the  pressure  balance  at  the  boundary  layer  between  plasma  and  neutraJ  gas 
in  a  gas  divertor  of  tokamak  reactors. 
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5  A  Two-Optical-Path  Laser  Fluorescence  Method 

Measuring  the  Ha  laser  induced  resonance  fluorescence  scattering  from  plasma  is  a  method  to 
determine  the  neutral  h\drogen  density.  This  metht>d  has  been  applied  to  toroidal  tokomak 
plasma  il7  ,  18],  19i.  and  tandem  mirror  plasma  ^20  .  21!. 

The  laser  induced  fluorescence  signal  is  usually  much  weaker  then  the  Ha  emission  from 
plasma  and  the  laser  stray  light.  The  signal  extraction  is  difficult  and  inaccurate.  This  work 
shows  that  a  clean  signal  can  be  extracted  by  using  a  two-optical-path  technique.  One  f)ptical 
path  was  focused  on  the  spot  in  the  plasma  shone  by  the  laser  light,  therefore  it  carries  all 
signals  (the  Ha  emission  from  plasma,  laser  stray  light  and  fluorescence  signal).  Another  path 
is  focused  on  the  spot  inside  the  plasma  slightly  off  the  laser  shined  spot,  It  therefore  carries 
only  the  Ha  signal  and  stray  light  which  will  be  called  background  signal  for  convenience. 
The  optical  paths  were  separated  by  a  prism  and  then  led  to  two  photomultipliers.  By 
subtracting  the  electric  signals  from  the  photomultipliers  using  a  differential  amplifier  we 
were  able  to  eliminate  the  background,  and  a  nearly  pure  fluorescence  signal  was  obtained. 


5.1  Experimental  Method 

This  Two-Optical-Path  (TOP)  laser  fluorescence  diagnostic  system  was  located  in  the  south 
exhaust  chamber  of  the  machine  at  a  distance  of  25  cm  from  the  closest  mirror  point  of  the 
end  cell.  The  electron  temperature  and  density  were  measured  by  a  fast  Langmuir  probe 
22\ 

The  crossectional  view  of  the  \'acuum  chamber,  the  laser  beam  and  the  detection  optical 
system  is  shown  in  Fig.  21.  The  laser  used  here  is  a  flashlamp-pumped  dye  laser  (Candela 
EDL-6)  tuned  by  an  angle-tuned  40/x  air-spaced  etalon.  A  monochromator  was  used  to 
monitor  the  wavelength  of  the  laser  output.  A  photo-diode,  mounted  directly  behind  the 
99.9%  reflecting  re«r  cavity  mirror  of  the  laser  is  used  tr)  monitor  the  power  of  each  laser 
pulse.  A  telescopic  system  in  front  of  laser  improves  the  beam  convergence  .  The  beam  can 
be  focused  at  various  radial  position  in  the  plasma  by  adjusting  the  concave  lens.  In  the 
beam  entry  port  there  were  three  copper  knife-edge  baffles  w'hich  were  arranged  to  minimize 
the  stray  light  intensity.  The  gate  valve  in  front  of  the  baffles  was  used  for  vacuum  isolation. 
The  scattering  system  housed  in  a  box  on  top  of  \acuum  chamber  can  be  moved  radially. 
As  will  be  discussed  later  it  contained  a  big  collection  lens,  7  cm  in  diameter  and  10  cm 
focal  length,  mounted  across  a  30cm  x  10cm  observation  window.  Radial  scan  can  be  made 
from  7-  =  0  to  r  =  4.5  cm.  This  limitation  was  due  to  the  10  cm  transverse  opening  of  the 
aperture  which  can  be  increased  for  larger  radial  scan  in  the  future.  A  30cm  x  10cm  viewing 
dump  which  was  made  of  stacked  non-magnetic  razor  material  located  on  the  bottom  of  the 
chamber  opposite  to  the  scattering  system  box.  The  beam  dump  at  the  exit  port  consists 
of  two  disks  of  absorbing  OB  10-type  glass.  The  glass  was  set  at  Brewster's  angle  to  the 
incident  l2iser  direction  in  such  a  w'ay  that  both  polarization  components  were  attenuated. 

The  sketch  of  the  scattering  system,  the  light  paths  and  side  view  of  the  plasma  and 
vacuum  chamber  are  shown  in  Figure  25.  The  laser  is  normal  to  the  entry  quartz  window 
and  shines  on  the  plasma  at  the  spot  as  shown.  The  optical  paths  consist  of  a  big  lens,  a 
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scattehinc; 


Figure  2-1;  The  laser  beam  and  deleelion  opliral  svhtein. 


TO  DATA  SYSTEM 


Figure  2.5:  The  sketch  of  the  scattering  system 
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filter,  a  prism  and  two  small  lenses.  The  surface  of  the  prism  is  deposited  with  a  thin  coating 
of  aluminum  which  gives  a  surface  reflectivity  of  The  prism  was  adjusted  so  that  the 

image  of  the  each  optical  path  fell  on  the  entrance  slits  of  the  photomultipliers,  PM  I  and 
PM2  respectively.  PMl  receives  the  light  from  the  laser  shone  plasma  spot,  whereas  PM2 
receives  the  light  from  the  plasma  spot  at  an  axial  position  slightly  off  from  the  laser  spot. 
To  induce  the  fluorescence  emission  of  //«  from  the  hydrogen  the  wavelength  of  the  laser  was 
tuned  to  6563.4.  The  induced  fluorescence  light  emitted  form  the  plasma  at  the  laser  spot  is 
led  to  PMl.  whereas  PM2  picks  up  the  background  signal  only.  When  the  laser  light  is  not 
in  tune  with  the  Ha  line  no  fluorescence  emission  was  induced  and  both  photomultipliers 
PMl  and  PM2  receive  only  the  background  signal.  The  outputs  from  P.Ml  and  PM2  are 
nearly  equal  and  can  be  subtracted  completely.  The  difference  of  the  signals  was  obtained 
using  a  differential  amplifier  (Tektronix  502)  with  a  gain  of  10.  The  contribution  of  the 
background  was  further  reduced  by  placing  a  narrow  band  filter,  with  a  bandpass  of  5  A  and 
transmission  coefficient  of  53%,  in  between  the  big  lens  and  the  prism.  The  Hamamatsu  R928 
photomultiplier  tubes  is  shielded  magnetically  with  high  material.  The  optical  parts  were 
also  shielded  both  optically  and  magnetically.  The  output  from  the  differential  amplifier  and 
the  power  monitoring  signal  were  processed  by  6  MHz  and  100  kHz  digitizers  respectively. 

The  laser  was  tuned  to  Balmer-alpha  line  6563  A  by  an  etalon  and  had  a  linewidth  of 
about  0.75  .A.  The  laser  output  power  was  about  80  kW  and  1/3  of  that  power  was  focused 
«m  the  plasma.  The  focused  spot  size  was  0.4  — 0.5/cm^.  The  laser  power  of  10//jVTcm"^.\ 
was  sufficient  to  saturate  the  n  =  2  to  n  =  3  transition  of  hydrogen  atom.  The  fluorescence 
intensity  does  not  depend  on  the  laser  power. 

The  laser  fluorescence  signal  was  absolutely  calibrated  by  Rayleigh  scattering  in  nitrogen 
gas.  The  ratio  of  the  fluorescence  signal,  F,  to  Reyleigh  scattered  signal.  R,  is 

F  ^ 

where  A.V3  is  the  difference  in  population  of  hydrogen  atom  in  the  n  =  3  state  with 
and  without  laser,  A32  —  4.4  x  lO^sec'^  is  the  induced  transition  probability  from  n  =  3 
to  n  —  2.  cTn  —  2.16  X  lO'^^cm^  is  the  Reyleigh  scattering  cross  section,  Tf,  =  5ps  is  the 
duration  of  the  laser  pulse,  r  =  0.5  cm  is  the  radius  of  laser  spot,  =  3.54  x  10**F  (  P 
is  the  pressure  of  nitrogen  in  Torr)  is  the  molecular  number  density  of  nitrogen  in  Reyleigh 
scattering,  Nl  —  7.9  x  10'^  is  the  input  photon  number  for  Reyleigh  scattering  laser  light 
of  energy  Ei  =  480  mJ  used  in  our  experiment.  The  Reyleigh  scattered  signal  was  found 
to  be  R  =  7.8F.  Thus,  AiVa  =  4.49  x  10'’F,  where  F  is  fluorescence  signal  related  to  the 
Reyleigh  scattered  signal.  According  to  the  theory  of  the  collisional-radiative  modeI[9],  the 
ground  state  population  of  neutral  hydrogen  depends  on  the  electron  temperature  and  plasma 
density.  They  were  measured  by  a  fast  Langiniur  probe  located  at  same  radial  position  as 
the  laser  spot. 
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5.2  Results  and  Discussion 


The  flat  top  of  the  ICKII  pulse  duration  is  20.  The  laser  was  fired  at  12  ms  after  the  initiation 
of  K'KH  pulse  and  was  at  about  the  middle  of  flat  top.  Figure  26a  shows  the  light  spectrum 
from  path  P.\ll.  The  sharp  line  contains  the  laser  stray  light  and  fluorescence  light,  and  the 
continuous  spectrum  is  the  //„  emission  and  stray  lights.  Figure  26b  shows  the  spectrum  of 
the  lights  from  the  path  PM2.  The  difference  of  the  magnitudes  of  the  two  sharp  lines  is  the 
fltiorescence  signal  as  shown  in  Figure  27.  The  magnitude  of  this  real  signal  is  much  smaller 
than  the  stray  laser  light  and  even  smaller  than  the  emission  and  stray  lights.  Without 
the  use  of  PM2  the  extraction  would  be  very  difficult  and  inaccurate.  .As  is  shown  in  figure 
27  the  TOP  method  gives  a  clear  signal  free  of  noise. 

The  radial  profile  of  electron  temperature  and  density  measured  by  Langmuir  probe 
are  shown  in  Figures  28a  and  28b.  The  neutral  hydrogen  density  profile  .Vo,  inferred  from 
A.Va  with  electron  temperature  T,  and  density  is  shown  in  Figure  28(c).  The  electron 
temperature  T^is  about  40-50eV’  and  density  is  2  —  3  x  10® /cm^  which  is  an  order  of  mag¬ 
nitude  lower  than  that  in  the  center  cell  because  the  plasma  was  expanded  in  the  exhaust 
and  the  mirror  ratio  was  10.  The  neutral  density  near  the  center  is  one  order  of  magnitudes 
smaller  than  that  at  r  —  4.5.  The  plasma  is  fully  ionized  on  the  axis. 


56 


Figure  26:  The  signal  through  path  PMl  (Ha  ,  background  light  and  stray  laser  light  and 
the  signal  through  path  PM2  (  The  laser  induced  fluorescence  light  the  stray  light,  //„  and 
background  light. 
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Figure  27:  Fluorescence  light  spectrum  produced  by  the  TOP  method. 
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6  Experimental  Study 

The  experimental  setup  is  shown  in  figure  2.  The  wave  is  launched  from  the  antenna  in  the 
high  field  region  and  propagates  toward  the  resonance  plane  and  inward  toward  the  axis. 
I'he  measurements  were  made  at  the  resonance  plane  where  the  diagnostics  are  currently 
located. 

6.1  Diagnostics 

In  order  to  successfully  carry  out  the  experiments  described  in  here  we  have  to  continuously 
develop  new  diagnostics.  Following  are  the  ones  that  were  developed  this  year: 

1 .  \  retarding  field  grid  energ\-  analyzer  has  been  built  and  tested  successfully.  How'ever, 
funding  limitations  have  prevented  us  from  acquiring  the  necessary  function  generator 
and  power  supply  to  bring  this  unit  on  line. 

2.  A  wave  propagation  B-dot  probe  array  consisting  of  five  evenly  spaced  detectors  has 
been  built  and  installed.  The  axial  profile  of  the  B-field  measurement  was  made  using 
this  array. 

3.  A  triple  Langmuir  probe  for  electron  density  and  temperature  measurements  has  been 
built,  tested  and  has  been  used  for  experiments.  A  careful  calibration  with  single  probe 
measurements  has  to  be  made. 

4.  An  entirely  new  E-pr<»lie  has  been  successfully  developed.  This  makes  the  measurement 
of  the  electric  field  component  of  the  wave  possible  for  the  first  time. 

5.  A  laser  fluorescence  system  to  measure  neutral  density  has  been  mauie  operational  and 
has  been  used  for  the  exhaust  study. 

6.  A  Doppler  broadening  spectral  meter  has  been  developed.  This  system  wull  be  used 
for  ion  temperature  measurement. 

6.2  ICRF  Wave  Experiment 

The  experimental  setup  is  shown  in  figure  29  (Note  the  slightl  different  illustration  from 
Figure  2).  The  wave  is  launched  from  the  antenna  in  the  high  field  region  as  indicated  and 
propagates  toward  the  resonance  plane  and  inward  toward  the  axis.  The  measurements  were 
made  at  the  resonance  plane  where  the  diagnostics  are  currently  located.  The  magnetic 
component  of  the  wave  is  measured  by  a  B-dot  probe  array.  The  electric  component  is 
measured  by  an  E  probe.  The  density  is  measured  by  both  Langmuir  probe  and  microwave 
interferometer.  The  numderical  simulation  region  is  shown  by  the  shaded  area. 

The  electron  density  of  the  plasma  at  two  different  ICRF  frequencies  was  measured  with 
a  microw'ave  interferometer  and  is  shown  in  figure  30.  The  characteristics  of  the  density 
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results  shows  the  phenomena  of  the  existence  of  two  tvpes  of  arges  observed  from  the 

emission  6-  Figijre  30a  shows  the  discharge  for  a;  >  u-v,  and  is  classified  as  type  1. 
whereas  figure  30b  shows  the  discharge  for  u;  <  and  is  classified  as  type  II.  Type  II 
discharge  is  more  quiescent  than  type  I. 

Type  I  discharges  generate  densities  that  are  lower  than  type  II  discharges.  The  peak 
line  density  obtained  from  the  interferometer  is  approximately  3.2  ±  0.85  x  10**cm“*,  which 
corresponds  to  a  plasma  density  of  approximately  1.6  *  0.42  x  10”cm~^.  The  average  line 
density  is  approximately  2.6  ±  0.85  x  i0**cm~^  or  1.3  :r  0.42  x  lO'^cm"^.  The  uncertainty  is 
due  to  both  pickup  as  well  as  the  resolution  of  the  phase  digitizer. 

The  radial  profile  of  of  the  wave  for  the  the  type  II  discharge  measured  by  a  single 
B-dot  probe  at  midplane  is  shown  in  figure  31.  The  solid  curve  is  the  numerical  results  and 
the  open  circles  with  crosses  are  the  experimental  data.  They  are  in  good  agreement. 

The  axial  profile  of  of  the  ICRH  wave  for  the  type  II  discharge  measured  by  the  B-dot 
probe  array  is  shown  in  figure  32a.  It  is  seen  that  the  wave  is  damped  as  it  travels  toward 
the  resonance  plane.  The  numerical  simulation  of  the  wave  was  presented  previously  [6]  and 
a  trimetric  view  of  the  calculated  B,  is  shown  in  32b.  The  corresponding  region  of  the  axial 
profile  measured  by  the  B-dot  probe  array  is  indicated  in  figure  32b  as  the  damped  region 
between  the  peak  and  resonance  plane.  The  wave  is  rising  to  a  peak  and  is  then  damped 
toward  the  resonance  plane  in  agreement  with  the  experimental  data.  Such  a  damping  effect 
is  also  predicted  by  the  analytic  solution  discussed  in  the  theoretical  analysis  section. 

Tt*  determine  a  range  of  stable  operating  conditions,  a  series  of  experiments  at  various 
magnetic  fields  was  performed  to  see  the  effect  of  the  magnetic  field  intensity  on  the  stability 
of  the  plasma.  Instability  occurs  at  low  field.  Figure  33  shows  triple  Langmuir  probe  data 
illustrating  this  effect.  From  top  to  bottom,  the  plasma  stability  improves  when  the  field 
intensity  increases  (as  is  indicated  by  the  increasing  values  of  coil  current). 

6.3  Propulsion  Parameters 

It  is  our  plan  to  carry  out  a  series  of  of  experiments  by  making  simultaneous  measurements 
with  all  the  diagnostics  on  line.  To  do  so  the  machine  has  to  be  modified  to  increase  the 
number  of  ports  and  more  data  acquisition  equipment  have  to  be  acquired.  In  this  experiment 
we  try  to  utilize  our  existing  resources  to  obtain  the  badly  needed  propulsion  parameters.  In 
order  to  extract  the  propulsion  parameters  we  need  to  know,  at  least,  the  electron  2uid  ion 
temperature,  the  electron  and  ion  density,  the  radiation  losses,  the  gas  pressure,  the  mass 
flow-rate  and  thrust.  These  measurements  were  taken  simultaneously  with  a  fast  Langmuir 
probe,  a  triple  Langmuir  probe,  a  Doppler  broadening  spectral  meter,  radiation  monitors, 
a  residue  gas  analyzer,  fast  pressure  gages  and  the  ballistic  thrust  balance  just  developed 
for  this  purpose.  A  typical  plasma  shot  is  shown  in  Figure  34.  The  top  two  signals  are 
triple  probe  measurements  and  the  third  signal  is  electron  saturation  current  of  a  Langmuir 
probe.  The  Ha  radiation  from  hydrogen  atoms  is  shown  in  Figure  34  (fourth  signal)  without 
scanning.  The  bottom  signal  in  this  figure  is  the  Doppler  broadened  line  scanned  by  the 
rotating  mirror.  This  line  was  obtained  with  entrance  slit  of  the  spectrometer  set  wide  open, 
in  order  to  admit  enough  light  to  see  the  line,  thereby  the  resolution  is  very  poor.  The 
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Figure  29:  Simplified  experimental  setup  for  ICRF  wave  propagation  and  heating  experi¬ 
ment.  The  shaded  area  is  the  numerical  simulation  region 
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Figure  30:  Plasma  density  measured  by  interferometer  for  (a)type  I  discharge  (u;  >  Wci),  and 
(b)  type  II  discharge  (u;  <  Wci) 
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Figure  31.  Radial  profile  for  type  II  discharge 
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Figure  33:  Sleability  study 
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flediirerl  hydrogen  temperature  is  not  reliable.  In  order  to  obtain  a  high  solution  spectral 
line  we  have  to  deliberately  i  rc  ate  a  railiative  pla>ma.  Also,  for  the  purpose  of  measuring  the 
ion  temperature  we  have  to  seed  the  hydrogen  plasma  with  helium.  .Again  the  radiation  from 
the  helium  ions  is  very  weak.  We  have  not  yet  found  a  means  to  enhance  the  radiation  strong 
enough  to  permit  a  meaningful  measurement.  This  is  a  strong  indication  of  low  radiation 
losses  from  the  plasma.  The  Plasma  and  Propulsion  parameters  are  listed  in  Table  2  and 
their  implication  has  been  discussed  in  the  progress  summary  section. 
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Figure  3-1:  Sample  probe  signals  of  the  plasma  discharge:  (a)  triple  pn*l)e  voltage  signal, 
(b)  triple  probe  current  signal,  (c)  the  Langmuir  probe  saturation  c  urrent  signal,  the  fourth 
signal  (d)  and  fifth  signal  (e)  are  the  7/„  radiation  and  |)op|)ler  broadened  line  sjMctrum. 


7  Flight  System  Study 

We  have  undertaken  a  flight  system  design  study  and  mission  analyses  (postulating  a  piloted 
Mars  mission  at  10  MW’e  and  50%  efficiency)  11  .  1  he  preliminary  results  are  very  attrac¬ 
tive.  Figure  35  shows  the  extrapolated  operational  envelope  for  this  system.  Igp  \'ariations 
between  2,000  and  20,000  s  can  be  obtained  while  the  available  thrust  varies  between  10.000 
N  and  10  for  hydrogen  propellant.  This  is  a  unique  feature  for  this  propulsion  system. 
The  thrust  is  very  large  at  low  I,p.  and  decreases  as  the  vehicle  gains  speed  (low  thrust  at 
high  Isp)-  Such  a  capability  can  not  be  found  in  most  other  types  of  propulsion  schemes 
where  high  thrust  is  always  associated  with  high  I^p-  Therefore  the  vehicle  can  be  flown  at 
constant  acceleration  for  continuous  optimization  of  propulsive  efficiency.  The  thrust  can 
be  further  increased  by  using  heavy  elements  such  as  argon  or  cesium  as  propellant  at  low 
temperature.  The  major  advantage  of  cesium  is  that  it  can  be  preionized  with  waste  heat 
from  the  p«>wer  generator  and  other  equipment  such  as  rf  transmitters. 

Figure  36  illustrates  one  possible  system  layout.  For  a  10  MW  flight  system  the  rocket 
itself  has  a  specific  weight  of  0.04  kg/kW  while  the  rf  system  is  4  kg/kW.  Using  an  exist¬ 
ing  space  nuclear  electric  system  design  from  Rocketdyne  [23]  we  arrive  at  an  overall  figure 
of  8.04  kg/kW  making  the  variable  rocket  very  competitive.  The  general  subsystem 
specific  weight  breakdown  is  shown  in  Table  3.  These  preliminary  calculations,  based  on 
reproducible  experimental  results,  open  the  way  for  a  more  in-depth  subsystem  design.  This 
effort  is  ongoing  albeit  limited  by  manpower  constraints. 

Initial  calculations  performed  recently  for  our  group  by  J.  Gilland  and  C.  Hack  of 
NAS.A/Lewis  [I2]  indicate  an  attractive  150  day  transit  to  MARS.  The  total  trip  time  is 
480  days  including  a  90  day  stay  on  MARS.  The  calculations  used  an  input  v-alue  of  8 
kg/kW  and  an  efficiency  of  50%.  Since  then,  experimental  results  and  hardware  design 
modifications  have  allowed  us  to  refine  the  subsystem  specific  weight  estimates,  tis  well  zis  to 
consider  a  more  realistic  efficiency  \'alue  of  70%. 

From  the  mission  analysis  results  shown  in  Figure  37,  one  sees  the  variation  during 
both  the  accelerating  and  decelerating  phases  of  the  flight  leading  to  the  capture  of  the  ve¬ 
hicle  in  the  MARS  orbit  in  all  cases.  Thrust  system  design  refinement  using  state-of-the-art 
buperconductiug  technology  as  well  as  a  more  careful  evaluation  of  the  structural  require¬ 
ments  in  a  space-borne  system,  have  led  to  a  much  more  compact  and  lightw’eight  device. 
This  evolution  is  shown  in  Figure  38  which  compares  the  experimental  device  to  the  opti¬ 
mized  flying  system.  The  experimental  device  is  made  from  copper  and  stainless  steel.  The 
coils  are  discrete  and  large,  therefore  heavy  support  structure  is  necessary,  whereas  for  the 
flying  system  the  coils  are  sriperconductor  and  the  chamber  and  coil  casing  are  aluminum. 
The  overall  size  of  the  coils  is  reduced  by  a  factor  of  four.  Therefore  the  system  is  much 
more  compact  and  weighs  only  one  fifth  of  our  experimental  device.  A  conceptual  picture  of 
the  rocket  is  shown  in  Figure  39. 
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Tabic  3.  Spcrifii  Weight  «)f  The  Subsystems 

Plasma  Rocket  0.01  kg/kWe 

RF  Transmitter  1.0  kg/kWe 

i  -Nuclear  Electric  P^mer  4.0  kg/kWe 

pOverall  specific  weight  8.04  kg/  kWe 
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EXPECTED  PERFORMANCE  ENVELOPE 
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Figure  33;  An  extrapolated  operational  envelope  for  a  10  MWe  propulsion  system. 
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Pinj  =  6  XIW,  Hydrogen 
Efficiency  =  (iOfo 
EeX'haust  radius  =  5  cn 


Varying  Specific  Impulse  Optimization 
One-way  Mars  Mission,  10  MWe,  2016 
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(»)  Experimental  Device. 

Copper  magnet,  S.S  chamber. 
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Figure  38:  A  10  MW  'I'arulein  Mirror  I'lasina  Propulsion  SysUMii:  (a)  F.xpfrinu-ula 


8  Exhaust  Modeling 

To  model  the  exhaust  a  time  dependent,  three  fluid  (electron,  ion  and  neutral)  code  by 
solving  transport  equations  in  cylindrical  get»metry  including  interactions  between  species 
has  been  developed.  We  began  to  learn  the  dynamic  behavior  of  the  plasma-gas  system. 
The  theory  and  the  simulation  method  are  described  below. 

The  macroscopic  quantities  (n,V,  and  T)  are  modeled  via  a  fluid  representation.  The 
basic  equations  can  be  obtained  from  the  work  of  Braginskii  24l.  Braginskii  assumed  that 
there  was  no  processes  that  convert  particles  of  one  species  to  another,  so  the  right  hand  side 
of  his  particle  conservation  equation,  for  example,  is  set  to  zero.  We  have  such  processes,  and 
must  generalize  the  equations  appropriately.  We  need  fluid  equations  for  the  conservation 
of  particles,  momentum,  and  energy.  The  equations  to  be  presented  will  be  for  the  plasma 
ion  species.  The  equations  for  electron  species  are  identical.  The  equations  for  the  neutral 
gas  are  the  same,  except  for  the  momentum  conservation  equation,  which  lacks  the  Lorentz 
force  term.  To  close  the  fluid  equations,  we  need  an  equation  of  state  to  specify  the  relation 
between  pressure  and  temperature. 


8.1  The  Transport  Equations 

The  equation  for  conservation  of  particles  is  (Fq.  (1.11)  of  Braginskii  on  page  225  in  reference 

'241) 

^  V  •  (nV)  =  S,  (33) 

Here  n  is  the  particle  density  and  V  is  the  mean  fluid  velocity,  defined  as 

V  =  (v)  (34) 

V  is  the  exact  particle  velocity,  and  the  ( • )  denotes  an  average  over  velocity  space.  Sp  is 
the  particle  source  term,  to  be  described  in  detail  later. 

The  equation  for  the  conservation  of  momentum  is  (Eq.  (1.12)  of  Braginskii  in  reference 

n 

^  (mnV)  +  V  •  (mn  (vv))  -  en  -i-  -V  x  R  (35) 

Here  v  =  v'  +  V,  the  exact  particle  velocity,  and  v'  is  a  random  velocity. 

Ignoring  viscosity,  this  becomes: 

(mnV)  +  V  •  (mnW)  4-  Vp  —  en  -i — V  x  b')  =  R  (36) 

where  m  is  the  mass  of  a  particle,  c  is  the  speed  of  light,  p  is  pressure,  E  is  the  electric  field 
and  B  is  the  magnetic  field.  R  is  the  source  of  momentum,  to  be  described  in  detail  later. 
Note  that  both  of  these  equations  are  in  conservative  form,  and  this  is  the  way  that  they  are 
numerically  differenced  in  the  code. 
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For  conservation  of  energ\-.  we  start  from  the  heat-balance  equation,  (Eq.  (1.23)  of  Bra- 
ginskii  in  reference  21I) 

i  v  (^nTv)  +n7V  V  t  (37) 

Here  q  is  the  heat  flux  and  11  is  the  viscosity  tensor.  Q  is  the  heat  source  term,  to  be 
described  in  detail  later.  Ignoring  viscosity,  this  becomes 

2^^  T  V  •  (^nrv)  +  nrv  •  V  ^  V  •  q  -  g  (38) 

The  second  term  can  be  divided  as  follows: 


TV  •  nV  -  nV  •  VT  (39) 

The  first  term  above  can  be  combined  with  the  appropriate  time  derivative  term  of  Eq.  (8) 
to  give  the  continuity  equation  scaled  by  T. 

Writing 

Q-=Q-ItS,  (41) 

and  dividing  Eq.  (8)  by  (3/2)  n  leaves  us 

^  ^V.vr+  ^TV.V=  -(Q*- V.g)  (42) 

The  second  and  third  terms  can  be  rewritten  as 

v.(rv)-^rv-v  (43) 

o 

This  leads  to  the  following  equation  for  temperature: 

^+V.(rv)-irv-v=^(e--v-q)  (44) 

Heat  conduction  is  dominated  by  the  electrons,  and  thus  q  is  only  present  in  the  electron 
temperature  equation,  or,  q  =  q*.  The  form  for  qe  is: 

q.  =  KlVx(fcre)  (45) 


Here  k  is  Boltzman’s  constant,  k'  is  the  perpendicular  electron  thermal  conductivity  and  is 
given  as 


k'  =4.7 


nkTe 


(46) 


Here  u)ee  is  the  electron  gj-rofrequency  and  r*  is  the  basic  collisional  time. 
To  close  the  fluid  equations,  the  equation  of  state  for  pressure  is 

p  =  nkT 


(47) 
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8.2  The  Source  Terms 


The  source  terms  are  from  vari<»us  collision  processes  with  the  neutral  gas,  including  atomic 
reactions.  The  details  of  the  various  source  terms  will  now  be  presented. 

The  atomic  reactions  included  here  are  ionization,  charge-exchange  and  recombination 


and  are  characterized  as  follows. 

Electron  Impact  Ionization: 

f  X  ^  2c 

(48) 

Resonant  charge  -exchange: 

hot  cold  hot 

H*  ^  H‘‘ 

(49) 

Collisional-radiative  recombination: 

e.  //’  — »  H’  hi/ 

(50) 

and  electron  three-body  recombination: 

e~H^+e^H  +  e 

(51) 

In  the  recombination  reactions  above,  the  names  describe  the  way  the  energy  released  by 
the  recombination  is  carried  away.  H'  is  an  excited  state,  and  hi/  represents  photons. 

The  particle  source  term,  Sp.  consists  of  sources  from  atomic  reactions  between  the  plasma 
and  the  neutral  gas,  which  can  be  broken  into  three  terms 


9-9- 


i- 


S  -  5 


(52) 


Here  5ion  is  the  source  of  particles  from  electron  impact  ionization,  5cx  is  the  source  of  parti¬ 
cles  from  resonant  charge-exchange  and  is  the  source  of  particles  from  the  recombination 
reactions. 

The  specific  particle  source  terms  have  the  following  form: 


5rx  =  n,n2(«7v)rx  (53) 

nj  and  n2  are  determined  by  the  species  involved  in  the  reaction.  {(rv)„  is  the  reaction-rate 
for  the  specific  reaction.  The  reaction-rates  have  been  previously  calculated  and  stored  in  a 
lookup  table  to  be  used  by  the  code.  The  cross-section,  <r,  is  temperature  dependent.  The 
average  reiiction-rates  {<rv)  was  calculated  sissuming  a  Mzucwellian  velocity  distribution  for 
both  species  involved  in  the  reaction. 

The  momentum  source  term  has  the  same  form  as  the  particle  source  term.  Specifically, 

R  —  /?ion  +  Rex  ~  Rnc  (54) 

The  R's  represent  the  source  of  momentum  from  the  various  atomic  reactions,or: 

Rrx  =  (mV)j  n,n2(o’i’)rx 


1 1 


(55) 


fn  aHrlitifin,  the  R  term  mntains  elastic  collisions  between  the  plasma  and  neutrals,  i.e.. 
collisions  that  do  not  lead  to  atomic  reactiims.  These  elastic  collisions  are  commonly  referred 
to  as  friction.  The  frictional  momentum  transfer  can  be  modeled  as  (for  ions  scattering  off 
of  neutrals) 

i;)t'm  (56) 

where  is  the  collision  frequency  for  scattering  of  ions  by  neutrals. 

From  momentum  conservation  we  know  that  =  /2f,ni-  Electrons  also  scatter  elasti¬ 
cally  from  the  neutrals,  but  in  the  hybrid  model,  this  is  taken  into  account  in  the  term  for 
the  plasma's  resistivity  (to  be  described  later). 

The  last  source  term  to  consider  is  the  source  term  for  internal  energy,  Q.  This  includes 
the  internal  energy  brought  in  via  atomic  reactions, 


Q  —  Q'ion  Qcx  Qrct: 

(57) 

Each  Q  has  the  form 

<?rx  =  (^nkT^  ^nin2{cTv)r^ 

(58) 

For  the  ions  there  is  an 

additional  term  for  ion  heating  from  the  electrons 

3n2>A  nfc 

(59) 

Qi  =  -^-(Te-T) 

fUi  Te 


Electromagnetic  Terms 

All  that  is  left  is  to  find  equations  for  E  and  B.  vln  a  hybrid  plasma  model,  both  come 
friun  the  electron  m<»mentum  equation. 

In  the  limit  of  small  electron  inertia,  the  electron  momentum  equations  reduced  to 

VnkTf  1  „  , 

E  = - Ue  X  B  -i-  r/J  (60) 

en  c 

where  r)  is  the  plasma  resistivity.  The  last  term  of  this  equation,  tjJ,  is  simply  the  R  term  in 
the  notation  used  for  the  ion  momentum  equation.  In  other  words,  the  friction  of  electrons 
scattering  off  of  ions.  As  mentioned  above,  this  term  must  also  account  for  the  scattering  of 
electron  off  of  neutrals.  This  is  taken  into  account  in  the  rj  term. 

Vtolal  ~  Vtn  Vei  (hi) 

The  electron  momentum  equation  can  be  considered  the  Ohm’s  Law  for  the  system.  It 
should  be  emphasized  that  the  velocity  in  the  Lorentz  term  of  this  equation  is  the  electron 
velocity,  Ue-  not  the  fluid  velocity  V.  These  two  velocities  are  related  through  the  current, 


J  =  en(V  -  Ue) 

(62) 

Ue  =  -  ^  J  +  V 

(63) 

en 
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Thus,  the  II,  'X  B  term  lan  be  expanded  to  the  form 


V  X  B  ^  J  X  B  (6d) 

cn 

This  second  term  is  commonly  referred  to  as  the  Hall  current  term.  In  a  more  familiar  form, 
the  electron  momentum  equation  becomes 

E  +  V  X  B  =  qj  -  --  (J  X  B  -  Vp,)  (65) 

cn 

However,  we  will  use  Eq.  60  and  Maxwell's  equations  to  transform  the  electron  momentum 
equation  into  an  equation  for  B. 

The  radiation-free  limit  of  .Ampere’s  law  is 


47r 

V  X  B  =  —  J 

(66) 

c 

1 

(67) 

Faradav’s  law  is 


Introducing  the  magnetic  vector  potential  ( V  x  A  =  B)  in  the  Coulomb  gauge,  and  utilizing 
Maxwell's  equations  we  can  transform  the  electron  momentum  equation  to 


A  =  r^  ^UeXB--^77(VxB)  (68) 

en  47r 

The  9  component  of  this  equation  will  give  us  bf»th  B,  and  Finally,  if  we  take  the  curl 
of  the  electron  momentum  equation,  and  again  utilizing  .Jaxwell  s  equations,  we  obtain  the 
equation  which  was  used  to  advance  B$  numerically. 

The  electron  velocity  components  are  given  via  the  relationship  between  the  fluid  velocity 
and  current.  Specifically, 


Uer  - 
Itgr  — 
Uez  = 


1 _ 

'  4iren  dz 

+  r^V^A 

4iTen 


19, 
1  airB.) 


(38) 


8.3  Boundary  Conditions 

To  complete  the  formulation  of  the  problem  one  must  specify  consistent  boundary  conditions 
{BC's)  for  the  equations  discussed  above.  The  geometry  is  a  cylindrically  symmetric  duct. 
BC's  must  be  specified  at  the  inflow,  outflow,  r-  =  0,  and  the  wall,  r .  .o,.. 

The  inflow  is  assumed  completely  specified  for  the  plasma  and  held  fixed  in  time.  The 
same  is  assumed  for  the  neutral  gas  for  axial  injection.  If  radial  injection  is  used,  Zmin 
becomes  an  outflow  boundary  for  neutrals,  also  to  be  described  later. 
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At  the  origin,  r  0,  all  quantities  are  assumed  symmetric  in  6  which  leads  to  (he 
B('  di'dr  0.  This  holds  fur  all  quantities  except  4^.  The  BC  for  .4^  is  Dirichlet  zero. 
L'Hopital’s  rule  is  utilized  to  ensure  the  fields  derived  from  the  vector  potential  is  appropriate. 

At  the  wall,  r^ai,  4#  is  specified.  For  the  fluid  quantities,  most  are  set  as  Neumann  zero. 
The  exception  is  temperature.  In  some  cases  the  temperature  is  specified. 

The  most  difficult  BC's  occur  at  Smai-  Conventional  outflow  boundary  conditions  are 
used  for  the  neutral  gas,  e.  g.,  simple  interpolation  back  in  z.  \\  hen  the  same  BCs  are 
attempted  for  the  plasma,  the  computation  failed  eventually  due  to  the  problem  that  the 
plasma  fluid  reflects  (numerically)  from  the  outflow  boundary.  .Many  different  BC's  were 
attempted  for  the  plasma  fluid  at  2^0*-  The  best  results  were  obtained  by  simply  interpo¬ 
lating  back  up  a  flux  surface  (magnetic  field  line),  and  setting  the  fluid  quantities  constant 
along  the  field  line  just  at  Zmax-  In  this  sense  the  magnetic  field  line  is  approximating  a  char¬ 
acteristic,  even  though  in  actuality  it  is  not.  To  date,  the  best  results  have  been  obtained 
using  the  above  BC's  for  the  system. 


8.4  Discussion 

Our  problem  is  many  degrees  more  general  than  a  simple  one  fluid  MHD  problem  and  is 
therefore  many  degrees  more  complex.  Instead  of  solving  one  set  of  fluid  equations,  we 
are  solving  three  sets  of  fluid  equations  -  one  set  for  ions  and  one  set  for  electrons  plus 
the  neutral  transport  equation.  Instead  of  steady  state  we  are  solving  the  time  dependent 
problem.  The  reactions  between  species  were  taken  into  account.  The  coordinate  system  is 
cylindrical  instead  of  a  two  dimensional  Cau-tesian.  The  electron  stream  and  the  ion  stream 
arc  governed  by  different  time  scales.  The  plasma  time  scale  is  many  orders  of  magnitude 
faster  than  the  neutral  jet. 

The  expanded  exhaust  section  of  the  mirror  device  is  shown  in  Figure  40  and  is  marked 
by  the  coordinate  axis  r  —  2,  where  r  =  0  is  the  centerline  of  the  device.  Z  =  0  is  the  entrance 
of  the  plasma  into  the  exhaust.  The  flowing  plasma  paterns  for  various  gas  injection  method 
are  presented  in  Figure  41.  Only  lower  half  (below  the  z-axis)  of  exhaust  is  displayed.  The 
plasma  is  seen  flowing  from  left  to  right.  Figure  41a  is  the  flow  of  plasma  without  neutrals 
at  3.75/is.  The  plasma  is  seen  fanned  out  along  the  field  lines.  Figures  41b  show  that  plasma 
flow  with  coaxial  neuteral  jet.  The  plasma  was  contained  and  boundary  layer  was  formed. 
However,  the  neutral  jet  was  eroded  downstream  and  the  plasma  eventually  spread  out  along 
the  field  lines  (This  effect  is  not  shown  here).  The  reason  i  >bvious  because  the  neutrals 
can  not  flow  as  fast  as  plasma.  A  steady  state  flow  can  not  be  maintained.  Figure  41c  shows 
density  contour  of  plasma  and  neutrals  when  the  neutral  was  injected  radially  as  indicated 
by  dashed  lines  and  an  arrow.  The  plasma  was  contained  better  but  neutrals  were  eroded  at 
the  tip.  Figure  41d  shows  the  radial  injection  at  increasing  neutral  density  and/or  jet  width. 

In  the  radial  injection  case  the  flow  characteristics  w’as  further  improved  by  increeising 
the  width  of  the  radial  gas  jet.  The  flow  time  has  been  extended  to  23/i  which  is  nearly 
steady  state. 

These  analyses  show  that  injection  angle,  gas  flow  rate  and  density  are  very  important 
and  need  to  be  studied  thoroughly.  There  is  a  wfide  range  of  parameter  space  that  has  to  be 
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(a) 


i 

(b) 


Figure  -iO:  Plasma  exhaust 
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ill  orrl<T  to  understand  the  fiyriami«-  pi'iperties  i)f  this  olasma-gas  system.  However, 
tfie  plasma  iletai  hntenl  with  tht'  injei’tion  nf  annual  lold  gas  is  ch'arly  demonstrated.  The 
detathment  can  be  further  improved  with  a  tlisconnet ted  magnetic  system.  .\  disconnected 
magnetic  system  has  been  designed  and  the  detachment  has  been  demonstrated  with  single 
particl  model.  This  work  will  be  reported  when  rigorous  theoretical  analysis  being  done. 
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9  Future  Work 


To  conclude  this  research  it  is  useful  to  to  give  some  projections  and  outline  of  possible  fu¬ 
ture  plan.  To  do  so  some  remark  on  the  future  need  of  space  propulsion  and  the  capabilities 
of  this  system  is  in  order.  A  new  class  of  advanced  propulsion  devices  is  required  to  meet 
future  space  exploration  goals  including  the  efficient  delivery  of  large  payloads  to  geostation¬ 
ary  orbit  and  fast  manned  and  unmanned  missions  to  the  Moon  and  Mars.  The  variable 
l,p  tandem  mirror  rocket  provides  such  capabilities  by  operating  at  high  power  levels  (>  10 
MW)  and  high  efficiency  (>  60%).  The  system  is  also  capable  of  a  variable  thrust//,p  profile 
which  optimizes  the  vehicle  propulsive  efficiency  over  the  entire  flight. 

Other  distinctive  features  of  this  device  further  enhance  its  rele\'ance  as  a  viable  propul¬ 
sion  method.  For  example,  its  magnetically  insulated,  electrodeless  design  greatly  amelio¬ 
rates  materials  constraints  thus  enabling  high  power  operation.  Systems’  analyses  also  point 
to  a  favorable  specific  mass  scding  (  <  6kg/kW  at  10  MW)  for  a  full-up  superconducting 
design  coupled  to  a  nuclear  Rankine  power  conversion  system.  Additional  gains  are  envi¬ 
sioned  for  designs  operating  at  higher  voltages  (>  10  kV)  and  high  frequencies  (>  500  Hz). 
The  present  experimental  device  has  operated  virtually  trouble-free  for  over  four  years  and 
continues  to  serve  as  a  '‘pathfinder”  to  explore  new  design  approaches. 

In  spite  of  these  advances,  multiple  deficiencies  in  our  knowledge  base  continue  to  exist 
for  this  system.  For  example,  although  high  temperature/low  density  (high  7,p/low  thrust) 
operation  has  now  been  demonstrated,  it  remains  to  be  shown  that  this  system  can  run 
efficiently  in  the  high  density/low  temperature  (high  thrust/low  f,p)  regime.  Moreover, 
the  well-behaved  plasma  conditions  observed  thus  far  need  to  be  verified  when  the  device 
operates  asymmetrically  in  the  flowing  plasma  mode.  In  addition,  magnetically  vectored 
plasma  jets  eventually  have  to  detach  from  the  guiding  field  to  provide  useful  thrust.  While 
such  detachment  has  been  shown  successfully  in  numerical  models,  it  has  yet  to  be  confirmed 
experimentally.  Finally,  the  physics  of  ICRF  wave  propagation  in  tandem  mirror  plasmas 
has  not  been  fully  understood  despite  the  achievement  made  in  this  program,  such  knowledge 
is  crucial  to  maximizing  the  overall  power  absorption  efficiency  of  the  system  over  its  wide 
operating  range. 

As  mentioned  repetitively,  we  would  like  to  perform  the  experiment  to  demonstrate  the 
plasma  detaclimeiil  necessary  to  bring  this  tandem  mirror  rocket  to  the  mainstream  of  elec¬ 
tric  propulsion  program.  However,  resource  limitations  have  prevented  us  from  operating 
in  the  streaming-plasma  mode.  This  capability  requires  the  manufacturing  of  a  vacuum 
interface  to  the  large  vacuum  tank  in  ou.  laboratory  and  the  design  and  manufacture  of 
a  propellant  injector.  Also,  two  additional  magnet  power  supplies  are  required  to  run  the 
system  asymmetrically.  The  future  system  enhancements  are  highly  desirable.  The  near 
term  and  long  term  research  plans  are  presented  below\ 


84 


9.1  Near-Term  Research 

F«»r  the  near  future  we  should  continue  with  our  enveli»pe  expansion  appr<*ach  by  validating 
theory  «>ver  a  wide  operational  range.  Therefore  the  immediate  next  step  is  to  carry  out  the 
experiment  to  increase  thrust  by  increasing  the  density  and  reducing  the  temperature  and  to 
improve  the  energy  conversion  efficiency.  The  major  equipment  for  this  experiment  is  a  higher 
frequency  microwave  transmitter.  For  this  purpose  a  surplus  1-4.5  GHz  microwave  transmitter 
has  been  acquired  and  installed.  Therefore  no  major  capital  investment  is  required  except 
some  system  upgrade.  Plasma  parameters  in  the  central  cell  will  be  evaluated  not  just  at 
one  point  but  as  profile  functions  of  axial  and  radial  position.  These  data  will  continue  to  be 
used  to  generate  engineering  design  guidelines  applicable  to  spacecraft  systems.  In  particular, 
an  accurate  characterization  of  efficiency  and  its  dependence  on  plasma  conditions  will  be 
sought.  At  the  exhaust,  the  direct  measurement  of  thrust  and  1^  will  be  obtained  through 
two  detectors  presently  being  installed:  a  grid  energy  analyzer  to  measure  axial  ion  energy 
and  ballistic  balance  to  measure  total  thrust.  The  specific  tasks  are: 

1.  Continue  envelope  expansion  by  v'alidating  theory  and  experiment; 

2.  Obtain  profile  measurements  of  plasma  parameters  and  characterize  efficiency  and  its 
dependence  on  plasma  conditions: 

.3.  Obtain  direct  I^p  and  thrust  measurements  at  the  exhaust  with  grid  energ.v  analvzer 
and  ballistic  balance. 


9.2  Long-Term  Research 

It  is  fundamental  to  the  proof-of-principle  to  be  able  to  operate  this  device  in  the  asym¬ 
metric,  flowing-plasma  mode.  Such  an  upgrade  requires  two  major  hardware  modifications: 
construction  and  installation  of  a  vacuum  interface  to  the  vacuum  plenum  already  available 
in  our  laboratory  and  the  acquisition  of  two  additional  power  supplies  to  allow  independent 
magnetic  field  programming  at  the  two  end-cells.  Moreover,  with  completion  of  these  modi¬ 
fications,  we  would  also  propose  the  design  and  installation  of  the  plasma  injection  system, 
and  of  the  annular  nozzle  in  order  to  fully  test  plume  detachment  characteristics. 

In  the  theoretical  arena  we  plan  to  continue  our  numerical  code  expansion  and  validation, 
as  well  as  our  ongoing  efforts  on  mission  design  and  flight  system  analysis.  The  specific  goals 
are  as  follows: 

1.  To  achieve  the  asymmetric,  flowing-plasma  mode  of  operation; 

2.  To  work  toward  high  density,  low  temperature,  high  power,  high  thrust  operation. 

3.  To  make  four  major  heu'dware  modifications: 

•  Vacuum  interface  between  system  and  existing  plenum; 

•  Upgrade  of  rf  power  transmitters; 
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•  two  .ulditional  power  supplies  needed  for  independent  magnetic  field  programming 
at  the  en«l  «ells. 

•  r)esign  and  installation  of  fueling  injection  system  and  of  annular  nozzle  for  flowing 
plasma  and  Carry  <»iit  experimental  plume  detachment  studies. 

4.  Continue  the  ongoing  effort  on  code  expansion  and  validation  on  the  ICRF  wave  prop¬ 
agation  and  absorption  for  a  wide  range  of  plasma  parameters,  as  well  as  mission  and 
flight  system  design. 
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